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U.S.Army  Regional  Aerosol  Transport  Simulation 
(ARATS) 

Contract  DAJA45-90-C-0051:  Third  Interim  Report. 


Introduction. 


The'Sccond  Interim  Report  concluded  by  identifying  3'possi^'Iines  of  continuation: 

1.  The  influence  of  synoptic  type  upon  dispersion. 

2.  The  stability  of  multi-particle  models  used  for  long  range  transports,  and 

3.  Casework  amd  model  intercomparisons  associated  with  events  in  the  Gulyas  determined 

in  discussion  with  Dr  Holt  and  other  ARATS  scientists.  _ -  ^  '' s 

However^ the  adventitious  occurrence  of  the  Gulf  oil  fire  plumeTias  led  Di  Hol^to  redirect 
the  main  thrust  of  the  investigation  away  from  1  and  2  and  towards  3,  that  is,  studies  of 
the  spread  of  pollutants  from  Kuwait  using  a  range  of  numerical  modelling  techniques. 
At  the  2nd  ARATS  Workshop  it  was  agreed  between  Dr  Holt,  Dr  ApSimon  of  Imperial 
College  and  myself  representing  the  U.K.  Met  Office  thaHhe  continuation  and  conclusion 
of  the  work  under  the  existing  contracts  should  take  the'Torm  entirely  of  a  comparative 
study  of  the  numerical  modelling  of  the  plume  during  an  interval  when  it  was  particularly 
well  observed  by  aircraft  and  satellite,  using  Imperial  College’s  ‘3-DRAW’  model  and  the 
Met  Office  ‘NAME’  model.  These  are  both  long  range  transport  models  using  multiple 
particle  methods,  but  are  differently  structured  and  parametrized.  A  single  data  base  is  to 
be  provided  by  the  Met  Office.  Such  a  comparison  should  prove  invaluable  at  highlighting 
those  aspects  of  the  dispersion  parametrizations  which  can  lead  to  significant  differences 
in  simulation.  ^ _ 
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This  study  will  then,  at  Dr  Holt’s  request  (and  perhaps  also  in  recognition  of  the 
changing  geopolitical  scene),  to  replace  the  work  on  synoptic  type  over  Europe.  However, 
good  progress  had  already  been  made  on  the  2nd  of  the  above  projects,  a  detailed  study 
of  the  stability  of  the  NAME  model  with  respect  to  the  particle  numbers  used.  This  gave 
a  first  assessment  of  the  optimal  particle  numbers  required  by  a  ‘Monte  Carlo’  multiple- 
particle  model  adequately  to  simulate  the  spread  of  a  continental-scale  plume  over  a  period 
of  a  week  to  10  days.  This  may  well  be  the  first  time  this  has  been  attempted. 


Under  head  3  must  be  included  the  major  Met  Office  contribution  to  the  2nd  ARATS 
Workshop.  A  description  is  given  below  of  the  extensive  modelling  work  carried  out  in 
preparation  for  this  meeting,  supplemented  with  more  detail  in  Annexes,  and  of  the  Met 
Resesirch  Flight  contribution. 
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The  2nd  Workshop  on  the  Army  Regional  Aerosol  Transport 
Simulation- ARATS  2. 


Tlie  success  of  this  meeting  owed  much  to  the  efforts  of  U.K.  Met  Office  staff.  The 
following  attended: 


Dr  P  W  White,  Deputy  Director  (Short-range  Forecasting  Research). 

Dr  F  B  Smith,  senior  Met  Office  dispersion  scientist. 

R  H  Maryon,  Head,  of  Special  Investigations  and  Nuclear  Accident  Response. 
Mrs  S  P  Ballard,  meso-scale  NWP  modelling. 

D  W  Johnson,  Met  Research  Flight. 


On  the  first  morning  R  Maryon  and  Mrs  Ballard  presented  modelling  studies  of  the  Gulf 
oil  fire  plume  for  material  released  from  OOOOUCT  on  1st  March  1991:  a  case  chosen  in 
advance  by  Dr  Holt.  A  description  of  the  products  of  the  meso-scale  model  was  written  up 
by  Mrs  Ballard  and  R  T  H  Barnes  (Annex  1),  and  sample  output  from  the  corresponding 
integrations  of  the  NAME  regional  model  provided  by  R  Maryon  (Annex  2).  These  models 
are  complementary,  the  meso-scale  handling  the  near-source  plume  which  is  interacting 
with  the  local  meteorology,  NAME  the  wider  dispersion  problem.  NAME  utilised  smalyses 
of  the  meteorology  from  the  U.K.Met  Office  operational  numerical  weather  prediction 
(NWP)  models  while  the  meso-scale  model  was  only  initialised  therefrom,  and  preceded 
to  generate  its  own  developments  on  the  meso-scale. 

Both  the  models  simulated  the  plume  well  during  the  48  hr  integrations,  as  fax  as 
can  be  deduced  from  satellite  observations.  Both  picked  up  the  complex,  3-dimensional, 
sheared  structure  of  the  plume.  The  meso-scale  model  generated  interesting  meso-scale 
circulations,  induced  by  the  land/sea  juxtaposition  and  by  the  effects  of  the  absorption  of 
solar  radiation  by  the  dense  smoke.  Cells  of  ascending  air,  capped  by  pyrocumulus,  were 
forecast,  in  general  agreement  with  satelfite  and  aircraft  observations.  Such  effects  as  the 
temperature  drop  beneath  the  smoke  (due  to  its  shielding  effects),  and  general  levels  of 
concentration  were  reasonably  well  handled. 

The  NAME  model,  of  course,  is  a  regional-scale  model,  which  does  not  provide  the 
resolution  or  detailed  physics  of  the  meso-scale  model.  In  particular,  it  has  no  parametriza- 
tion  to  represent  the  interaction  of  the  plume  with  the  meteorology  of  the  atmospheric 
boundary  layer  (ABL) — the  usually  turbulent  region  close  to  the  surface.  Aircraft  and 
other  observations  revealed  that  the  plume  modified  the  ABL,  reducing  the  turbulent 
energy  available  at  low  levels,  and  accordingly  reducing  the  mixing  of  smoke  into  the 
ABL  until  the  plume  eventually  beesune  sufficiently  diffuse  to  allow  renewed  penetration 
of  solar  radiation  and  generation  of  turbulent  kinetic  energy  (TKE).  The  NAME  model 
diagnoses  an  ABL  depth,  but  always  assumes  that  the  ABL  is  well-mixed  by  TKE  in 
daytime  conditions,  and  accordingly  spread  material  through  the  ABL  from  the  start. 

Despite  this,  the  NAME  model  reproduced  the  3-dimensional  structure  and  spread 
of  the  plume  into  the  wider  field  very  well,  and  it  was  considered  that  the  model  pos- 
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sessed  value  in  (1)  providing  a  benchmark  umulation,  as  a  reliable  and.  tested  regional 
scale  model,  and  (2)  in  providing  a  Lagranjpan  (as  opposed  to  Eulerian)  approach  to  the 
problem  of  modelling  the  transport,  spread  and  deposition  of  pollutants  or  obscurants. 
Unlike  the  meso-scale  model  the  NAME  model  was  able  to  provide  predictions  of  the 
time-integrated  concentrations  in  air  (‘dosage’)  smd  most  importantly,  of  sedimentation 
of  pollutants  to  the  surface  and  of  their  washout  by  rain. 

The  lofting  of  smoke  by  the  absorption  of  solar  radiation  was  also  modelled  by  the 
NAME  model,  but  in  a  different  way  from  the  meso-scale  model.  As  NAME  is  a  La- 
grangian  multi-particle  model,  the  heating  of  particles  by  the  flux-divergence  of  solar 
radiation  had  to  be  computed,  and  the  individual  particles  then  lifted  buoyantly  to  their 
equilibrium  level.  Output  from  the  lofting  run  is  distinguished  in  Annex  2  by  the  let¬ 
ter  B.  It  will  be  noted  that  in  comparison  with  the  non-lofting  run  (labelled  A)  the  air 
concentrations  in  the  ABL  are  substantially  reduced,  and  that  in  the  upper  layer  (to  850 
mb)  somewhat  increased.  Dry  deposition  of  carbon  to  the  surface,  and  time-integrated 
air  concentrations  in  the  ABL  are  also  reduced.  The  results — the  lofting  generally  limited 
to  about  1  km,  except  for  insignificant  qusmtities  of  material — are  qualitatively  similar 
to  those  of  the  Eulerian  meso-scale  model,  and  indeed,  this  mutual  support  enabled  an 
early  prediction  in  Browning  et  al  (1991)  that  the  stratosphere  (and  hence  global  climate) 
would  not  be  threatened  by  lofted  smoke. 

A  future  study  might  well  look  at  the  nesting  of  the  meso-scale  and  regional  models,  the 
former  to  handle  the  interactive  plume  near  the  source,  the  latter  the  spread  of  material 
further  idield.  A  most  interesting  suggestion  is  the  initialising  of  the  regional  model  from 
the  output  of  a  prior  integration  of  the  meso-scale  model. 

On  the  first  afternoon  a  presentation  was  made  by  D  Johnson  of  the  Met  Research 
Flight  (MRF),  who  participated  in  the  MRF  C130  aircraft  traversals  of  the  plume  during 
March  1991.  Measurements  of  the  local  met  conditions,  radiation,  ur  and  plume  chem¬ 
istry,  cloud  physics  and  aerosol  were  made  using  the  comprehensive  instrumentation  of 
this  aircraft.  The  near-  and  far-held  structures  of  the  plume  were  also  elucidated:  a  fuller 
account  is  given  in  Johnson  et  al  (1991)  which,  with  Browning  et  al  (1991),  is  appended 
at  Annex  3  for  convenience.  D  Johnson  also  provided  the  Working  Group  with  a  com¬ 
prehensive  and  very  useful  summary  of  the  flights  (Gibbs  et  al  1991),  not  reproduced 
here. 

With  the  exception  of  D  Johnson  all  members  of  the  Met  Office  party  brought  their 
varied  experience  and  expertise  to  the  discussion  and  recommendations  which  are  set  out 
in  the  conference  report. 


Stability  of  the  multiple-particle  model. 

The  stability  of  multi- particle  models  used  for  simulating  the  long  range  transport  and 
dispersion  of  pollutants  and  obscurants  (item  2  above)  was  addressed  during  the  summer 
of  1991,  after  the  heavy  workload  generated  by  the  Gulf  environmental  crisis  had  receded. 
This  study  may  have  been  unique  in  the  scale  of  the  integrations  carried  out  on  the  NAME 
model  to  establish  adequate  ‘validation’  against  which  practicable,  run-of-the-day  model 
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configurations  could  be  assessed.  For  instance,  integrations  were  made  in  which  6000 
particles /hr  were  released  and  carried  by  the  model  every  hour  for  S4O  hours,  and  their 
transport  and  diffusive  spread  simulated.  In  another  case  12000  particles/hr  were  released 
for  72  consecutive  hours,  and  the  simulation  continued  r.p  to  240  hr.  Integrations  such  as 
these  required  many  Mbytes  of  core,  and  the  runs  consumed  up  to  9  hr  of  CPU  time  on 
the  Met  Office  Hitachi  computer.  No  attempt  has  been  made  to  cost  these  overheads! 

Additional  code  was  written  and  added  to  the  model  to  carry  out  the  extensive  sta¬ 
tistical  analysis  required,  and  the  results  were  plotted  by  writing  graphics  routines  for  a 
Microvax  computer  and  transferring  the  data  between  the  systems.  The  results  are  unex¬ 
citing,  but  of  very  considerable  practical  importance:  it  is  valuable  to  have  some  idea  of 
the  optimal  number  of  particles  which  should  be  released  in  a  multi-particle  model  (a)  to 
simulate  a  given  situation  as  well  as  the  model  is  capable  of  doing  but  (b)  not  to  tie  up 
computer  storage  and  CPU  time  processing  large  numbers  of  particles  which  are  adding 
very  little  to  the  accuracy  of  the  simulation.  The  results  showed  that  for  a  Chernobyl- 
scale  incident  there  is  not  a  great  deal  to  be  gained  over  the  European  domain  from  re¬ 
leasing  more  than  SOO  particles/hr,  and  extremely  little  improvement  obtainable  beyond 
1000  particles/hr,  when  limitations  in  the  model’s  diffusion  parametrization  probably  put 
constriunts  upon  the  achievable  accuracy. 

Opportunity  was  taken  to  write  up  these  results  for  the  NATO/CCMS  19th  Interna¬ 
tional  Technical  Meeting  on  Air  Pollution  Modelling  and  its  Application,  held  at  lerapetra, 
Crete,  Sept  1991,  in  view  of  the  desirability  of  communicating  these  results  to  the  wider 
modelling  fraternity.  A  copy  of  this  paper  (Maryon  1992)  is  at  Annex  4. 


Plans  for  remainder  of  Contract. 

In  discussions  with  Dr  Holt  and  Dr  ApSimon,  Imperial  College,  it  was  agreed  to  complete 
the  project  with  a  detailed  and  structured  comparison  of  Imperial  College’s  3- DRAW  and 
the  Met  Office’s  NAME  long  range  transport  and  dispersion  models.  The  integrations 
would  start  at  OOOOUCT  27th  March  1991 — a  period  for  which  both  satellite  and  aircraft 
observations  of  the  Gulf  plume  are  available.  The  models  will  use  identical  release  profiles 
in  time  and  space,  and  in  source  emission  strength:  this  is,  for  some  reason,  an  ideal  that 
has  not  quite  been  attained  hitherto.  The  output  will  consist  of  trajectory  endpoints,  ABL 
and  elevated  level  concentrations  of  carbon  in  air,  and  wet  and  dry  deposition  of  carbon 
to  the  surface.  Data  for  this  analysis  is  to  be  assembled  and  provided  by  R  Maryon  at 
the  Met  Office.  Dr  ApSimon’s  and  her  colleague  Dr  Lowles’  more  detailed  proposals  ^ixe 
at  Annex  5.  As  discussed  in  the  Introduction,  a  comparison  of  two  similar  but  diff  .rently 
parametrized  models  using  identical  scenario  and  data  should  be  particularly  i  ivaluable 
in  highlighting  areas  of  potential  uncertainty  and  difficulty  in  this  type  of  simulation. 

The  relevance  of  this  work  to  ARATS  lies  not  only  in  the  assessment  of  the  suitability 
of  models  of  this  type  for  application  in  practical  situations,  and  the  reliance  which  cim 
be  placed  upon  their  products,  but  also  in  the  insights  which  can  be  gained  into  areas 
where  the  models  have  the  potentiality  of  improvement.  It  is  hoped  that  the  results  of 
this  study  will  be  avidlable  by  the  end  of  March  next. 
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ANNEXES: 


Annex  1:  An  account  of  the  application  of  the  Met  Office  meso-scale  model  to  the  Gulf 
plume,  commencing  OOOOUCT  Ist  March  1991,  with  examples  of  the  model  output. 
(Mrs  S  P  BaUard). 

Annex  2;  NAME  output  fields  of  trajectory  end-points,  concentrations  of  carbon  in  air, 
wet  and  dry  deposition  of  carbon  to  the  surface  and  time-integrated  concentration 
(dosage)  for  the  situation  of  Annex  1.  The  charts  labelled  A  are  without  lofting  of 
the  plume  due  to  the  absorption  of  solar  radiation;  the  charts  labelled  B  include 
this  effect. 

Annex  3:  Copies  of  the  1991  Nature  articles  by  Browning  et  al  and  Johnson  et  al. 

Annex  4:  A  copy  of  Maryon(1992). 

Annex  5:  Agreed  proposals  for  the  last  phase  of  the  work  under  contract,  by  H  M 
ApSimon,  Imperial  College,  London. 


R.H.  Maryon 

Project  Manager,  NAME. 

Met  0(P),  U.K.  Met  Office, 

London  Road,  Bracknell,  Berkshire,  RG12  2SZ. 
December,  1991. 
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1 .  Introduction 

Since  29th  January  1991  30  hour  mesoscale  model  forecasts  for  an  area  in 
the  Gulf  have  been  run  routinely  twice  dolly  off  special  fine-mesh 
forecasts  for  that  area  (  off  global  forecasts  since  4th  June  ).  As  the 
fine-mesh  model  grows  its  orography  in  the  first  3  hours  of  the  forecast  it 
was  sensible  to  start  the  mesoscale  forecasts  after  the  finemesh  fields 
were  in  balance  with  the  finer  orography  and  forecasts  are  started  at  T+6. 
Forecasts  are  run  with  purely  interpolated  fine-mesh  forecast  data  as  the 
initial  conditions  <  now  global  forecast  data).  For  these  forecasts  smoke 
sources  are  included  and  the  smoke  is  advected  as  a  passive  tracer  (  ie  the 
interaction  with  radiation  is  neglected  ).  For  special  investigations  the 
model  con  be  run  including  the  interaction  of  radiation  with  the  smoke. 

2.  Grid/fflodel  specifications. 

The  grid  was  originally  centred  on  lat/long  29. 5N  47. 5E  and  uses  a 
transverse  Mercator  projection  with  origin  at  centre  and  great  circle 
running  N-S  along  47. 5E. 

Grid  corners  were  NW  35.35N  40.23E  HE  35.35N  54.77E 
SW  23. SON  41.05E  SE  23.30N  53.95E 
No. of  grid  points  were  89  W-E  91  N-S  at  15km  spacing, 
i.e.  grid  covered  1320  x  1350  km. 

There  were  32  model  levels  (terrain  following)  at  1.25,  2.5,  5,  10,  20,  40, 
70,  115,  175,  255,  355,  480,  630,  810,  1020,  1265,  1550,  1880,  2260,  2695, 
3190,  3750,  4380,  5085,  5870,  6740,  7700,  8755,  9910,  11170,  12540,  14025m. 
Albedo  =  0.25  increasing  to  0.6  if  snow  covered  . 

From  25th  Februery  1991  to  4th  June  the  surface  resistance  to  evaporation 
has  been  set  to  500  sm-1  in  the  desert  areas  and  100  sm-1  over  Iran  and  the 
Euphrates  valley.  Previously  it  was  set  to  lOOsm-1  everywhere.  Since  4th 
June  interpolated  unified  model  values  have  been  used  outside  the  desert 
areas. 

Also  on  25th  February  it  was  decided  to  change  the  orography  from  the 
values  interpolated  from  the  fine-aiesh.  The  new  orography  is  interpolated 
from  a  U. 5. Navy  10'  resolution  global  orography  dataset  onto  the  mesoscale 
model  grid.  A  simple  interpolation  gives  maximum  increase  in  height  of 
880m  between  2  adjacent  grid  points  on  the  western  slope  of  the  Iranian 
mountains  and  a  max.  height  of  3077m,  and  this  much  steeper  and  higher 
orography  caused  one  case  with  strong  SWly  winds  to  blow-up  in  the  HE 
comer  of  the  nndel.  Therefore  it  was  decided  to  interpolete  in  such  a  way 
as  to  give  more  detail  over  low  ground  (below  lOOOm)  but  a  similar  level  of 
detail  to  the  flnc-mesh  over  the  mountains  of  Iran.  The  new  max.  height 
is  2638  m  and  the  max.  increase  between  adjacent  grid  points  440m.  All 
points  with  interpolated  orographic  heights  less  than  Om  are  set  to  -Im  and 
assumed  to  be  sea.  This  resolves  the  Gulf  well  without  the  need  to  modify 
by  hand. 

Roughness  length  is  wind  speed  dependent  over  water.  Since  25th  February 
new  values  have  been  used  over  land.  Previously  the  roughness  length  was 
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O.tm  for  heights  greater  than  SOOm  and  0.001m  elsewhere.  The  new  values 
over  land  have  been  Interpolated  from  1.25* long  x  0.833* let  data,  prepared 
for  the  unified  model,  onto  the  mesoscale  grid.  Low  roughness  length 
values  cover  the  desert  SW  half  of  the  area  with  little  dependence  on 
orographic  height  while  the  low  lying  but  cultivated  Euphrates  valley  has 
higher  values  than  before.  The  Iranian  mountains  have  high  roughness 
lengths  as  before. 

On  19th  April  the  area  was  moved  northwards  but  dimensions  kept  the  same. 
On  8th  May  the  area  was  expanded  to  cover  the  original  southern  area  as 
well  as  the  north  and  the  number  of  levels  were  reduced  to  21. 

The  model  formulation  is  as  described  in  Golding  1990  (  Met  Mag  119  pg  81- 
96).  Particular  attention  is  drawn  to  the  radiation  schemes  used  in  the 
model.  A  modified  form  of  the  longwave  radiation  scheme  of  Roach  and 
Slingo(1979)  is  included.  Heating  rates  ere  applied  every  timestep  at  every 
model  level.  The  heating  rates  are  derived  from  flux  divergences  which  are 
recalculated  every  15  minutes  .  The  transmission  functions  are  represented 
as  products  of  functions  for  water  vapour,  carbon  dioxide  and  cloud,  the 
latter  being  weighted  linearly  by  the  maximum  cloud  fraction  in  the  slab. 
The  short-wave  radiation  budget  is  recalculated  every  time-step  but  only  at 
two  levels,  the  ground,  and  the  highest  cloud  top  if  below  5km.  At  cloud 
top,  absorption  is  parameteri 2ed  as  a  function  of  solar  zenith  angle  only, 
and  is  constrained  not  to  exceed  the  long-wave  cooling  in  the  same  layer. 
At  the  ground  the  transmission  is  obtained  by  applying,  mutipl icatively , 
transmission  coefficients  for  Rayleigh  scattering,  water  vapour  absorption, 
aerosol  scattering  and  cloud. 

3.  Inclusion  of  smoke  plume. 

The  transport  of  the  smoke  is  modelled  by  using  the  equation 

6c/6t  =  -v.Vc  +  6(Khe>c/6n)/t>n  +  KVZc  +  Sc 

where  c  is  the  aerosol  concentration  in  kg/m/grid  square 

Kh  is  the  vertical  mixing  coefficient  for  heat  =  Ih  S'* 

=  fnf  1 ,E.N2 ,S2 ) 

K  =  22500  m2s-l 

Sc  =  source  in  kg/m/grid  squere/s 
all  boundary  conditions  are  c=0 
w'c'  for  the  surface  layer  =  -  Ch  lv|  c<l) 


dw'c’/6n  =  0  at  the  top  of  the  model 

Cain  =  5  X  10-7  kg/m/grid  square  is  lower  limit  on  concentration  to  remove 
any  negative  values  generated  by  the  finite  difference  scheme. 

Since  18th  April  the  concentrations  have  been  in  kg/kg  so  that  the 
transport  is  conservative.  However  there  is  little  difference  in  the 
resul ts. 


4.  Source  of  annke  pi  lane. 

In  the  routine  forecests  sources  of  smoke  were  originally  added  at 
locations  identifed  from  previous  satellite  imagery.  Depths  of  sources  are 
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also  estimated  from  the  imagery  and  the  preliminary  conclusions  of  the  Gulf 
Task  Group.  The  latter  estimated  an  initial  plume  rise  of  1 -3km  based  on 
typical  winter  Oanuary/February >  conditions.  This  is  consistent  with 
reports  from  the  Gulf  and  the  evidence  from  satellite  imagery  that  the 
majority  of  the  visible  smoke  is  remaining  within  the  boundary  layer. 
Therefore  for  most  runs  the  source  has  been  kept  below  3km.  Runs  with 
injection  up  to  5km  appeared  to  transport  too  much  smoke  with  the  upper 
level  winds.  Since  25th  February  a  depth  of  2260m  has  been  used  and  the 
forecast  smoke  locations  have  generally  agreed  well  with  the  imagery 
considering  the  limitations  of  the  method  of  specifying  the  depths  and 
locations  of  the  sources.  It  is  even  harder  to  specify  the  total  and 
relative  strengths  of  the  sources  .  This  is  why  the  routine  forecasts  are 
run  neglecting  the  interaction  with  radiation.  However  the  strengths  have 
been  estimated  based  on  the  imagery  and  using  an  upper  limit  of  5 
Megatonne/year  from  all  sources.  The  source  is  then  distributed  uniformly 
with  height  from  the  surface  to  source  top  over  a  single  grid  point. 

From  25th  February  to  20th  March  the  sources  were  kept  constant  as  3  point 
sources 


1 ocet i on 
30. 5N  A7.SE 
29. 3N  47. 6E 
29. 2N  47. aE 


depth  strength 

2260m  5/38  x  5Mtonne/year 

2260m  15/38  x  5Mtonne/year 

2260m  15/38  x  5Mtonne/year 


Better  locations  could  be  defined  using  the  locations  of  known  oil  fields. 
However  the  above  sources  provided  good  predictions  and  since  the  output 
was  provided  to  various  users  the  locations  were  not  updated.  Forecasts 
with  more  distributed  sources  provided  very  similar  results.  On  20th  March 
the  number  of  source  points  was  increased  to  J5  covering  the  main  oil 
fields,  each  point  had  3/38  x  5Mtonne/year .  This  gives  a  total  input  of 
greater  than  5Mtonne/year  but  oil  fields  in  northern  Kuweit  were  included 
and  it  was  not  certain  whether  they  were  included  in  the  initial  source 
estimate,  also  source  estimates  were  being  questioned. 


Although  forecasts  are  run  routinely  no  attempt  ,  except  in  some  special 
runs,  is  made  to  pass  smoke  from  forecast  to  forecast.  Therefore  it  is  not 
possible  to  compare  the  early  stages  of  the  forecast  with  satellite 
observations  as  the  smoke  starts  from  zero  in  each  run. 

The  plume  is  assumed  to  be  non-buoyant  and  smoke  removal  (  eg  by  washout 
and  dry  deposition  >  is  ignored  apart  from  that  implied  by  the  boundary 
conditions  and  surface  fluxes. 

5.  Inclusion  of  absorption  of  solar  radiation  by  smoke 


If  the  interaction  of  the  smoke  with  radiation  is  to  be  studied  the 
absorption  of  solar  radiation  by  the  smoke  can  be  modelled  as  suggested  by 
Slingo  (  personal  communication  )by  using 
SC2)=  So  cosCZ)  exp(-ku/  cosCZ)) 

where  2  is  the  solar  zenith  angle. 

Sc  cosCZ)  Wm- 2  is  *he  incoming  solar  radiation  at  the  top  of  the 
a  tmosphere 

k  =  absorption  coefficient  =  IC  m2  g- i 
n=32 

u  =  Z  C(n)  Ac 
n=n( z ) 
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CCn) 


concentration  at  leve 


n  in  gni-3  =  c(n)xl000/(  15x15x1x106  ) 


The  heating  rate  is  then  calculated  as 

66] /Ot  =  Cg/cpit>6S/6F' 

The  solar  fluxes  and  pressures  are  then  calculated  at  midlevels. 

In  the  implementation  of  the  scheme  the  absorption  due  to  smoke  is  used  to 
reduce  the  solar  fluxes  reaching  cloud  top  and  the  surface.  However,  in 
order  tc  simplify  the  inclusion  of  the  effects  of  smoke,  the  solar  fluxes 
used  in  the  calculation  of  heating  rates  due  to  smoke  absorption  are  not 
reduced  by  the  effects  of  clouds,  water  vapour,  scattering  and  background 
aerosol.  This  means  that  there  is  a  mismatch  between  the  surface  and 
atmospheric  heating  rates  but  the  smoke  absorption  is  a  worst  case  scenario 


6.  Heat  of  combustion 
This  is  not  included. 

7.  Forecast  for  02  3/3/91 


Charts  are  attached  f rom  the  30  hour  forecast  from  182  1/3/91  valid  at  02 
3/3/91.  These  are  from  the  routine  forecasts  with  passive  smoke  -  ie  no 
interaction  with  scler  radiation.  This,  forecast  has  not  been  studied  in 
detail  and  the  output  is  just  that  produced  routinely. 

Chart  A  shows  lOm  wind  and  present  weather  with  orographic  height  contours 
every  lOOOm.  f  =  visibility  less  than  1km,  -  =  visibility  l-Skm.  The 
visibility  does  not  include  the  effects  of  smoke  only  cloud  water,  humidity 
and  a  default  background  value  of  aerosol  suitable  for  the  UK  (  ie  not 
adjusted  for  desert  areas  etc  ) 


Chart  B  shows  l550m  wind  and  the  vertically  integrated  column  smoke  density 
in  kg/grid  square.  Horiiontel  lines  are  a  problem  with  the  contour  routine. 


Chert  C  shows  48Cm  wind  and  the  smoke  concentration  at  480m  in  kg/grid 
squere/m.  Horizontal  lines  are  a  problem  with  the  contour  routine. 


Chert  D  shows  1020m  wind  and  the  smoke  concentration  at  1020m  in  kg/grid 
square/m.  Horizontal  lines  are  a  problem  with  the  contour  routine. 
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Environmental  effects  from  burning 
oil  wells  in  Kuwait 
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Model  calculations,  constrained  by  satellite  obser¬ 
vations,  indicate  that  most  of  the  smoke  from  the 
oil  fires  in  Kuwait  will  remain  in  the  lowest  few 
kilometres  of  the  troposphere.  Beneath  the  plume 
there  is  a  severe  reduction  in  daylight,  and  a  day¬ 
time  temperature  drop  of  —10  X  within  —200  km 
of  the  source.  Episodic  events  of  acid  rain  and 
photochemical  smog  will  occur  within  —1,000- 
2.000  km  of  Kuwait.  But  changes  in  the  Asian 
summer  monsoon  are  unlikely  to  exceed  the  natural 
interannual  variability  and  stratospheric  ozone 
concentrations  are  unlikely  to  be  affected. 


The  main  environmental  effects  arising  from  the  burning  oil 
wells  in  Kuwait  come  from  emissions  of  oxides  of  carbon, 
sulphur  and  nitrogen,  and  unbumt  hydrocarbons  and  particu¬ 
lates.  The  scenario  used  for  this  study^  produced  aAer  dis¬ 
cussion  with  the  UK  Department  of  Energy  before  the  oil  wells 
had  been  set  alight,  assumes  80  Tg  of  oil  bunting  over  a  year, 
equivalent  to  the  annual  pre-invasion  production  of  about  1.5 
million  barrels  per  day,  with  6%  of  the  oil  being  converted  to 
fine  particulate  smoke.  Small^  assumed  1.6  million  barrels  per 
day,  with  8.6%  of  the  oil  being  converted  to  smoke.  The  products 
resulting  from  our  scenario  are  given  in  Table  1,  together  with 
some  current  national  and  global  emission  figures  for 
comparison. 

We  have  used  a  simple  model  of  plume  rise  to  investigate  the 
initial  behaviour  of  the  plume.  We  then  used  the  Meteorological 
Office’s  mesoscale  weather-prediction  model  to  estimate  the 
additional  increase  in  height  resulting  from  solar  heating,  the 
plume  shape,  and  the  daytime  reduction  in  surface  temperature, 
within  a  few  hundred  kilometres  of  Kuwait.  For  greater  dis¬ 
tances,  out  to  several  thousands  of  kilometres,  we  modelled  the 
plume  with  the  Meteorological  Office's  long-range  dispersion 
model,  which  also  provides  estimates  of  wet  and  dry  deposition. 
We  also  used  data  from  a  trajectory  model  as  a  basis  for  general 
circulation  model  integrations  to  determine  the  effect  of  the 
plume  on  the  Asian  summer  monsoon. 

Our  study  indicates  that  the  smoke  plume  is  very  unlikely  to 
reach  the  stratosphere,  and  Che  bulk  of  the  smoke  will  remain 
within  the  lower  troposphere  before  being  deposited  on  the 
ground  within  a  week  or  so  of  its  emission.  Consequently, 
changes  in  the  Asian  summer  monsoon  resulting  from  the  smoke 
are  likely  to  be  small  compared  with  the  natural  interannual 
variability.  Beneath  the  plume  within  about  200  km  of  the  source, 
it  is  likely  that  there  will  be  a  reduction  in  light  to  near  night-time 
levels,  and  a  reduction  of  the  daytime  maximum  temperature 
of  ‘*’iO*C.  The  plume  may  remain  visible  for  up  to  a  few 
thousand  kilometres,  and  slight  reductions  in  daylight  and  tem¬ 
perature  can  be  expected  beneath  it  at  distances  of  several 
hundred  kilometres.  Episodes  of  severe  add  rain  and  'black' 
snow,  and  photochemical  smog,  comparable  with  the  highest 


concentrations  encountered  around  major  industrial  areas,  are 
likely  to  occur  out  to  distances  of  a  thousand  kilometres  or  more 
downwind.  The  carbon  dioxide  release  (1%  of  current  global 
emissions)  will  have  an  insignificant  effect  on  global  warming. 

In  general,  we  agree  with  Small^  that  effects  will  be  significant 
in  the  Gulf  region  but  insignificant  on  a  global  scale. 

Local  effects 

During  the  first  few  minutes  after  emission,  materia)  in  a  smoke 
plume  will  rise  rapidly,  driven  mainly  by  the  heat  from  the 
burning  oil.  If  we  assume  that  the  initial  plume  rise  from  each 
well  is  not  enhanced  by  the  effect  of  other  wells,  and  use  a  heat 
source  based  on  that  expected  from  an  average  oil  well  (up  to 
500  MW),  then  simple  calculations  for  turbulent  plumes^  indi¬ 
cate  that  the  smoke  will  reach  an  initial  height  of  between  t 
and  2  km  in  the  range  of  stable  conditions  typical  of  winter. 
(Sma))^  suggested  a  height  of  700  m,  but  that  was  using  a  heat 
source  of  319  MW.)  In  summer,  when  the  well-mixed  surface 
boundary  layer  deepens,  smoke  can  be  expected  to  rise  to  the 
top  of  the  boundary  layer,  which  can  be  3  km  deep  (occasionally 
4  to  5  km).  In  unstable  situations,  that  is,  in  the  presence  of 
thunderstorms,  the  smoke  could  rise  significantly  higher,  and 
some  smoke  might  reach  the  upper  troposphere;  however,  analy¬ 
sis  of  the  Meteorological  Office's  climatological  archives  shows 
that  these  situations  occur  infrequently  and  are  generally  of 
short  duration. 

After  the  initial  rise  of  the  plumes  they  spread  out  and  merge 
together.  A  further  slow  rise,  termed  'self-lofting'^,  may  occur 
through  solar  heating.  The  oonsequence  of  this  effect  on  the 
overall  plume  was  investigated  by  using  the  Meteorological 
Office’s  mesoscale  numerical  weather-prediction  model^,  which 
has  a  gridscale  of  15  km,  adapted  to  include  transport  and 
diffusion  of  smoke  in  a  manner  similar  to  that  developed  for 
nuclear  winter  studies^  Transport  of  smoke  by  sub-gridscale 
convection,  and  removal  by  wet  and  dry  deposition,  were  not 
included.  In  this  and  the  other  numerical  models  described  later, 
incident  solar  radiation  was  assumed  to  be  absorbed  exponen¬ 
tially  by  the  smoke  with  an  absorption  coefficient^’*  of  10  m^  g~' . 
(Smalr  assumed  6.5m^g~'  with  a  slightly  greater  amount  of 
smoke.)  The  impacts  of  smoke  on  both  scattering  and  infrared 
radiation  were  ignored.  Nuclear  winter  studies^  show  that 
neglecting  the  infrared  properties  of  the  smoke  is  a  reasonable 
approximation  when  the  particle  sizes  are  predominantly  less 
than  one  micrometre.  The  small  infrared  effects  will  act  to 
ameliorate  the  predicted  daytime  surface  temperature  reduction 
and  may  reduce  nocturnal  surface  cooling.  The  reduction  of 
solar  radiation  by  overlying  cloud  was  also  ne^ected,  except  at 
the  surface,  to  simplify  the  model  and  to  investigate  the 
maximum  possible  impact. 

Initial  and  boundary  conditions  for  the  mesoscale  model  were 
interpolated  from  the  UK  Meteorological  Office's  operational 
regional  mode)  forecasts.  Throughout  the  period  of  the  forecast, 
smoke  was  injected  uniformly  from  the  surface  to  a  specified 
plume  top.  Forecasts  were  run  for  five  selected  days  in  January, 
February  and  March  1991  with  and  (as  a  control)  without  the 
interactive  effects  of  absorption  of  solar  radiation  by  smoke.  We 
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selected  cases  to  cover  the  two  main  synoptic  situations:  cloudy 
cyclonic,  when  the  smoke  was  advected  to  the  north  over  Iraq 
and  Iran,  and  clear  anticyclonic,  when  the  smoke  was  advected 
to  the  southeast  over  Saudi  Arabia  and  the  Gulf.  For  January 
and  February,  the  smoke  sources  were  distributed  uniformly 
over  a  square  of  side  45  km,  from  the  surface  to  a  height  of 
4  km.  For  the  later  cases,  point  sources  of  smoke  were  inserted 
from  the  surface  to  2.S  km,  at  locations  identified  from  satellite 
images.  The  plume  height  was  deduced  by  comparing  the 
forecast  vertical  wind  profile  with  the  direction  of  transport  of 
the  bulk  of  the  observed  smoke. 


PKk  1  Thirty-hour  forecast  of  the  Kuvraiti  smcM  pfume  vaHd  at  23 

February  1991  (see  also  Fig.  2c).  from  the  mesoscale  model.  «  Column- 
ifttegrated  dtotribution  of  smoke.  51  in  mg  m'*  (optical  depth-S  xir  xiO~^ 
where  the  absorptior)  coefficiem  k«lOm^g~'X  h  OffTerenoea  in  screen- 
level  (1.25  m)  temperature  due  to  absorption  of  ader  radtodon  (contours 
at  mtarvals  of  -5X  with  contour  addad  at  -2*C).  Tha  amai«  aouroa. 
centrad  at  29.1*N.  47.8* E.  is  5  magatons  par  year  spread  over  a  equare 
of  side  45km.  with  its  top  at  4km.  Tha  modal  domain  is  IJSOby  1.320km 
in  tha  horizontal  by  I4km  in  tha  vatHcaf  with  a  horfionm  resokition 
of  15km  and  32  lavals  spacad  irregularty  in  tha  vertical. 


TA8LE  1  Assumed  amuai  production  of  the  Kuvraiti  fires 

T^  of  emission 

Amount 
(U  pef  yr) 

Comparison  writh  current 
emtesion 

Fine  particulate  black 
smoke* 

5 

Roughly  a  tNrd  of  tire  carbon 
particles  produced  through  tropicat 
biomass  burning” 

Sulphur  (as  sugihur 
oxides) 

2 

Slightly  more  than  the  current  UK 
annual  sdphur  emissions” 

Nitrogen  (as  nitrogen 
oxides) 

0.5 

1968  (JK  nitrogwt  oxMe 
emissions”  were  0.75  Tg 

Carton  (ultimately  as 
carbon  dioxide) 

60 

About  1%  of  current  global  arwiual 
carton  dioxide  emissions  from 
fossil  fuel  combustion 

*  Small^  assumes  58  Tg  of  smoke  per  year, 
t  SmaK  (persor^  communication)  suggests  3.3%. 

The  ngures  are  compared  with  current  national  arwl  global  figures.  The 
percentage  of  sulphur  in  Kuwaiti  oil  vLtias  from  fleM  to  field,  but  an  average 
value  is  2.S%t  (figure  provided  by  the  Compartment  of  Energy).  The  percentage 
of  nitr^en  in  die  fuel  is  ^0.17%  (Dep^ment  of  Energy),  but  to  this  must 
be  added  the  quantity  of  nitrogen  oxides  produced  during  combustion  (above 
^2.000K);  we  have  assumed  an  overaH  figure  of  0.6%. 


In  the  control  runs,  the  smoke  remained  below  the  top  of  the 
initial  source  for  four  of  the  cases.  In  the  fifth  forecast,  south- 
vvesterly  flow  over  the  Iranian  mountains  lifted  significant  quan¬ 
tities  of  smoke  from  2.S  to  4  km.  and  trace  amounts  reached 
8-12  km  in  regions  of  convection.  Peak  concentrations  of  smoke 
were  forecast  just  above  the  top  of  the  boundary  layer  in  all 
forecasts  irrespective  of  the  initial  height  of  the  source. 

When  the  effects  of  solar  absorption  were  included,  lofting 
of  smoke  in  significant  quantities  was  limited  to  a  rise  of  1-3  km 
above  that  in  the  control  integrations.  With  flow  down  the  Gulf, 
and  solar  absorption  included,  small  amounts  of  smoke  reached 
8  km.  It  is  not  clear  whether  this  was  due  solely  to  selMofting 
or  partly  due  to  the  effects  of  induced  mesoscale  circulations, 
possibly  through  changing  land-sea  boundary  effects  (the 
model,  in  agreement  with  satellite  imagery,  sometimes  generated* 
cloud  in  the  regions  where  smoke  was  elevated). 

The  predicted  rise  of  the  smoke  is  low  compared  with  the 
typical  height  of  the  tropopause  (15  km)  over  Kuwait  Therefore, 
with  the  smoke  failing  to  reach  the  stratosphere,  its  lifetime  can 
be  expected  to  be  relatively  short  (see  later),  and  close  analogies 
with  nuclear  winter**"  and  large  volcanic  eruptions"  are  in¬ 
appropriate. 

We  also  assessed  the  impact  on  surface  temperatures  of 
absorption  of  solar  radiation.  There  is  a  large  variability  in  the 
position  of  the  plume,  and  the  associated  cooling  and  reduction 
in  light  levels,  because  of  the  day-to-day  variability  of  the  wind. 
Figure  lo  shows  the  smoke  distribution  (compare  also  with 
Fig.  2  later),  and  Fig.  1  b  gives  the  associated  reduction  of  surface 
temperature  for  23  February  1991.  The  maximum  daytime  tem¬ 
peratures  at  screen  level  are  depressed  by  about  10^  beneath 
the  hi|^  concentrations  of  smoke  within  about  200  km  of  the 
source,  and  temperature  drops  of  one  or  two  degrees  are  predic¬ 
ted  hundreds  of  kilometres  downstream  from  the  source.  The 
drop  of  20*C  shown  in  Fig.  \b  near  the  source  r^on  is  likely 
to  be  excessive.  The  model  caniMM  resolve  the  details  of  the  flow 
on  that  scale,  and  the  smoke  concentrations  may  be  large  enough 
that  infrared  effects  can  no  longer  be  ignored.  Results  will  also 
be  sensitive  to  the  assumed  surface  characteristics  and  the 
modelling  of  the  surface  layer  mixing.  Sunlight  is  reduced  to 
near  night-time  levels  in  the  model  when  thick  smoke  is  over¬ 
head,  but  the  model  cannot  represent  sub-gridscale  structure, 
and  this  may  allow  a  little  sunlight  to  penetrate. 

From  mid-March  to  mtd-Aprii,  we  have  compared  daily 
operational  forecasts  of  midday  temperatures  in  the  absence  of 
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the  smoke  with  such  observations  as  were  available  from  stations 
beneath  the  plume  (typically  4>8  observations  per  day  when  the 
plume  extended  south,  very  few  when  it  extended  north).  Checks 
on  the  forecast  values  were  provided  by  stations  in  unaffected 
regions.  This  subjective  assessment  suggests  that  within  *^250  km 
of  the  source  the  reductions  in  midday  temperature  under  the 
plume  ranged  from  5  to  8^,  decreasing  to  1-2  (probably 
the  limit  of  accuracy)  beyond  about  7S0km. 

Effects  at  longer  range 

The  Meteorological  Office's  long-range  dispersion  model'^  was 
used  to  investigate  the  longer-term  distribution  of  the  smoke, 
and  also  the  acid  rain  and  photochemical  smog  that  may  occur 
before  concentrations  reach  more  normal  background  levels. 
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The  model  simulates  the  spread  of  neutrally  buoyant  airborne 
pollutants  (for  example  fine  soot  or  sulphur)  by  releasing  very 
large  numbers  of  ‘particles’  into  the  model’s  atmosphere  at  the 
source  of  the  pollution.  The  particles  are  transported  by  the 
model  winds  and  vertical  motions,  and  are  diffused  vertically 
and  horizontally  by  random  perturbations,  which  represent  the 
small-scale  turbulent  motions  likely  to  be  experienced  in  the 
real  atmosphere  (especially  near  the  surface).  The  model  has 
been  run  daily  since  February  to  monitor  the  concentration  of 
the  smoke  particles  for  a  number  of  levels  in  the  vertical.  These 
Integrations  neglected  self-lofting. 

Six  integrations  were  re-run  with  solar  absorption  included, 
as  in  the  mesoscale  model,  and  the  particles  allowed  to  self-loft 
in  parcels  to  the  level  of  neutral  buoyancy.  Comparisons  between 
the  integrations  with  and  without  solar  absorption  showed  that 
in  strong  sunlight  a  small  proportion  (typically  5%)  of  the 
model’s  particles  lofted  rapidly  through  a  kilometre  or  more. 
Thinning  of  the  plume  by  entrainment,  coupled  with  the  model’s 
normal  turbulent  diffusion,  reduced  the  lofting  to  the  point  at 
which  other  vertical  motions  dominated,  usually  within  one  or 
two  days.  For  an  injection  height  of  2  km,  the  amount  of  smoke 
reaching  3-Skm  was  slightly  augmented,  but  most  of  the 
material  (90%)  remained  at  lower  levels. 

Close  to  the  source,  results  were  similar  to  the  mesoscale 
model.  An  example  of  the  shape  of  the  smoke  plume  at  different 
heights  is  shown  in  Fig.  2a  and  6.  In  the  boundary  layer,  below 
1.5  km,  the  plume  extends  southwards  and  fans  out;  between 
1.5  and  3  km,  it  encounters  winds  with  a  westerly  component, 
and  it  drifts  eastwards  across  the  Gulf.  Pictures  from  the  geo¬ 
stationary  satellite  Meteosat,  monitored  with  a  dedicated  com¬ 
puter  display  system,  broadly  confirmed  the  model  forecasts. 
The  visible  image  for  the  time  corresponding  to  Fig.  2a  and  6 
shows  the  plume  rapidly  broadening  towards  the  south  as  it 


no.  2  forty-oi^-hour  forecast  vaNd  at  12aMT.  23  February  1991.  from 
me  (ong-range  dtaperaior)  model,  ai  Shape  of  the  Kuwaiti  smoKe  plume  over 
the  hai^  Marval  0  to  13toyt  h  1.5  to  3  km.  aa  derired  from  the  long-range 
daperaion  modal  aaaiimlng  that  tha  amoha  axtandad  from  0.5  to  4km 
altituda  at  aourca.  Thaaa  may  ba  eomparad  wKh  tha  analyala  of  tha  Mereoaat 
vlatila  Imaga  for  tha  corranponifing  tkna  (Fig.  2c).  ebare  cross^wtehed. 
hatchad  and  abpplad  ahadng  repreaant  aublactivaly  dMarminad  ragkma  of 
dWfarent  optical  tWdmaaaaa.  Tha  rtadwrtframa  me  ahowa  tha  area  in  Fig.  1. 
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encounters  flows  at  different  levels,  before  thinning  to  such  an 
extent  as  to  be  undetectable  (Fig.  2c). 

To  estimate  the  lifetime  of  the  smoke,  the  ‘mass’  (correspond¬ 
ing  to  a  small  proportion  of  the  source  emission)  carried  by 
each  model  particle  was  progressively  reduced  to  take  account 
of  the  estimated  wet  and  dry  deposition  to  the  surface.  A 
standard  deposition  velocity  for  dry  particles  (O.OS  cm  s~',  valid 
over  the  boundary-layer  depth)  and  scavenging  coefficient  in 
rain  (].2x  ]0'‘*s~' permmh'')  were  assumed  for  die 
smoke'^'**.  Initial  and  boundary  conditions,  and  rainfall  rates, 
were  taken  from  the  operational  analyses  and  forecasts  of  the 
Meteofologicai  Office’s  regional  mode).  Fourteen  cases  were 
considered.  The  half-life  of  the  material  was  typically  S-IS  days, 
except  in  very  wet  conditions,  when  it  was  shorter,  and  in  deep 
dry  boundary  layers,  when  it  was  longer.  SmalH  estimated  a 
half-life  of  S  days  and  MacCracken'^  computed  3.5-6  days, 
depending  on  the  depth  of  the  release. 

Although  dry  deposition  is  the  principal  process  for  removing 
gaseous  sulphur  dioxide  from  the  plume,  it  is  less  important  for 
‘sulphur’  aerosol  or  sulphur  attached  to  particulates,  and  the 
time-integrated  dry  deposition  at  a  given  location  is  unlikely  to 
cause  significant  ecological  damage,  except  possibly  within  the 
immediate  area  of  the  fires  where  the  droplets  of  unbumt  oil 
fall  out  under  gravity.  Wet  deposition,  although  contributing 
less  to  the  overall  depletion  of  sulphur,  will  often  be  associated 
with  locally  very  acidic  rain.  Measurements  of  pH  in  eastern 
Saudi  Arabia '  ’  taken  before  the  Gulf  War  had  a  mean  of  almost 
5.5,  and  therefore  the  associated  background  ions  are  unlikely 
to  reduce  the  acidity  generated  by  the  plume  significantly.  We 
made  calculations  assuming  plume  widths  (typically  —200  km) 
and  depths  (3-4  km)  based  on  our  model  results  and  satellite 
observations.  Allowing  for  saturation  in  the  uptake  of  acidic 
species  into  growing  raindrops,  the  results  indicate  that,  once 
rain  is  encountered,  its  acidity  will  depend  only  weakly  on 
downwind  distance  out  to  at  least  2,000  km.  with  pH  values  in 
the  range  3.0-3.6  on  many  occasions.  These  very  acidic  values 
could  damage  crops  and  vegetation'^. 

Over  the  higher  mountains  much  of  the  precipitation  falls  as 
snow  and.  where  it  is  affected  by  the  emissions,  the  snow  may 
be  blackened  enough  by  soot  to  darken  it  considerably.  Forecasts 
from  the  mesoscale  model  (such  as  in  Fig.  In)  occasionally  show 
average  smoke-column  densities  of  -*~40mgm~^  over  the 
southern  Iranian  mountains.  If  this  were  scavenged  by  snowfall 
with  a  water  equivalent  depth  of  1  mro,  the  mass  mixing  ratio 


R6.  3  Mean  tftambutton  of  amoM  for  the  period  Aprt-Septente  dirtved 
from  lo-o^r  treiectories  comumed  using  the  Pe)Kter>  feoMy  of  the 
Europeen  Oentis  #of  Mitfum.nmge  tWmthei  Ttif imtlng**  The  ffgiw  rfwwi 
oonioure  of  the  number  of  hoera.  during  die  2-yser  period  Uffd-fO.  when 
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normtfiisd  by  die  umieapondiig  number  of  heurt  over  me  KuwsR 


of  soot  within  the  snow  would  be  typically  10~^,  which  could 
double  the  absorption  of  solar  radiation'  Spring  snow-meit 
could  therefore  occur  earlier  and  faster,  with  the  possible  con¬ 
sequence  of  high  run-off  and  highly  acidic  pulses  in  the  rivers. 

Close  to  the  source  of  the  plume,  the  low  light  levels  are 
expected  to  prevent  any  significant  photochemical  effects,  except 
at  the  top  of  the  plume.  As  the  plume  disperses,  however,  light 
levels  will  recover  to  values  more  typical  of  the  region.  This  will 
allow  rapid  oxidation  of  hydrocarl^ns  and  other  gases  (such 
as  HjS)  within  the  plume,  probably  over  a  timescale  of  a  few 
days,  resulting  in  locally  severe  episodes  of  photochemical  smog. 

Effect  on  the  monsoon 

To  determine  whether  the  smoke  emissions  are  likely  to  affect 
the  Asian  monsoon,  we  carried  out  an  idealized  sensitivity 
experiment  with  a  version  of  the  Meteorological  Office's  atmos¬ 
pheric  general  circulation  model  that  has  been  used  extensively 
for  seasonal  predictions  of  tropical  rainfall'^;  this  version  was 
chosen  because  it  produces  a  good  simulation  of  the  Asian 
monsoon.  The  radiative  effects  of  the  smoke  were  included  as 
in  nuclear  winter  studies'*;  this  approach  was  similar  to  the 
simple  formulation  in  the  mesoscale  model  described  above, 
except  that  we  included  the  effect  of  shading  by  clouds  on  the 
radiative  heating  of  the  smoke. 

The  distribution  of  smoke  used  in  this  model  was  derived 
from  the  wind  fields  actually  observed  in  two  previous  years, 
1989  and  1990,  using  the  global  trajectory  facility  attached  to 
the  numerical  weather-prediction  model  of  the  European  Centre 
for  Medium-Range  Weather  Forecasting'*.  A  smoke  particle 
was  released  once  a  month  from  each  of  the  vertices  of  an  area 
of  l*x  1.5”  covering  Kuwait  at  each  of  the  four  levels  950,  850, 
700  and  500  mbar,  and  their  three-dimensioDa)  trajectories  were 
followed  during  an  assumed  10-day  lifetime.  Figure  3a  and  b 
shows  the  resulting  density  of  smoke  at  two  different  levels, 
averaged  over  the  period  of  the  monsoon,  April  to  September. 

For  the  sensitivity  experiment,  we  prescribed  purely  absorbing 
smoke  (absorption  coefficient  10  g''),  with  a  smoke  loading 

of  0.012  gm'^  (vertical  optical  depth  0.12)  centred  on  Kuwait 
and  mixed  evenly  throughout  the  troposphere  up  to  a  height  of 
8  km.  The  distribution  of  the  smoke  in  Fig.  3  was  approximated 
by  aV(a^+  r^),  where  r  is  the  distance  from  source,  and  a  was 
set  at  two  model  grid  lengths  ( —600  km).  The  total  smoke  loading 
is  equivalent  to  assuming  a  smoke  lifetime  of  10  days  with 
emissions  at  a  rate  of  5  Tg  per  year. 


squwa.  for  tralectortM  a  balow  800  nmur  wid  6  between  800 «id 400  mbar . 
Tbe  centoura  may  be  roi#ey  tmerpretid  aa  normabrntf  conoantrabons. 
WthDu^  (hay  uniggrata  the  amount  of  amolit  at  long  rartgt  baoeuat  (hay 
do  not  ttha  mm  aooount  tot  tea  by  dapoaWon.  Tha  gap  and  laappaarmca 
of  maiartal  at  tongar  rants  am  due  to  ffia  ffiraa-dtoranatonal  nature  of  the 
tfstoemnaa  and  9m  imafl  sampm  tmed 


The  experiment  and  the  control  (with  no  smoke)  were  run 
from  April  to  September,  (n  the  experiment,  atmospheric  solar 
heating  in  the  centre  of  the  plume  increased  by  about  80  W  m  ^ 
(or  0.8 ’’C  per  day),  whereas  the  surface  heating  was  reduced 
by  *-25  W  m  V  The  surface  and  troposphere  remained  convec* 
tively  coupled,  so  there  was  net  heating  of  the  troposphere- 
surface  system  which  produced  a  slight  enhancement  of  the 
Asian  summer  monsoon  and,  in  most  places,  of  the  associated 
precipitation.  This  is  in  contrast  to  nuclear  winter  simulations 
in  which  the  smoke  is  so  thick  and  high  that  solar  heating  of 
the  surface  is  effectively  eliminated,  the  surface  and  troposphere 
are  thermally  decoupled^,  and  there  is  a  weakening  of  the  Asian 
summer  monsoon  with  reduced  precipitation^'. 
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Finally,  we  estimated  the  interannual  variability  of  precipita¬ 
tion  in  the  monsoon  region,  averaged  over  July  to  September, 
using  data  from  existing  simulations  run  with  the  observed  sea 
surface  temperatures  for  seven  different  years  ( 1950,  1958, 1976, 
1983, 1984. 1987  and  1988,  encompassing  wet  and  dry  conditions 
in  the  Sahel The  predicted  changes  in  precipitation  over 
southeast  Asia  (5-30*’  N,  70-105^  E)  exceed  twice  the  standard 
deviation  of  the  interannual  variations  over  only  10%  of  the 
area,  and  are  all  increases.  In  summary,  the  assumed  scenario 
does  not  lead  to  significant  decreases  in  the  precipitation  in  the 
Asian  summer  monsoon,  and  might  even  lead  to  increases  in 
certain  regions.  □ 
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Airborne  measurements  in  the  densest  part  of  the 
smoke  plume  at  about  120  km  from  the  burning 
wells  in  Kuwait  in  late  March  1991  showed  typical 
particulate  mass  densities  of  500-1.000  jig 
mixing  ratios  of  500-1.000  p  p.b.v.  of  sulphur 
dioxide  and  30-60  p  p.b.v.  of  nitrogen  oxides.  One 
thousand  kilometres  from  Kuwait,  ozone  con¬ 
centrations  in  the  plume  exceeded  background 
levels  by  about  50  p.p.b.v.  The  oil  burn  rate  was 
estimated  from  sulphur  fluxes  to  be  3.9  =  16  mil¬ 
lion  barrels  per  day  Significant  amounts  of  smoke 
were  observed  only  below  5.000  m  altitude,  and 
the  measured  attenuation  of  solar  radiation  by  the 
smoke  was  similar  to  those  assumed  in  recent 
assessments. 


ABOUT  600  naturally  pressurized  oil  wells  were  set  alight  m 
Kuwait  in  late  February  injecting  massive  ouamities  of 

smoke,  sulphur  dioxide  'SO;».  unburnt  hydrocarbons  and 
nitrogen  oxides  i  NO .  »  NO  NO: '  atmosphere  Assess¬ 

ments'"  of  the  impact  of  these  pollutants  are  subject  to  large 
uncertainties  in  the  magnitude  and  the  characteristics  of  the 
emissions  and  the  resulting  plume.  To  investigate  the  chemical 
and  particulate  composition  of  the  smoke,  and  its  effect  on 
visible  and  thermal  radiation,  the  C  l,^0  aircraft  of  the  Vfeteoro- 
logical  Office  Research  Flight  made  57  h  of  obsenations  in  eight 
flights  into  and  around  the  plume  in  March  19dl  Here  we  reoorx 


the  first  airborne  measurements  of  the  plume  both  close  to  the 
source  I  up  to  200  km » and  in  the  far  field  1 1.0<J0  km  from  source i 

Meteorology  and  plume  behaviour 

Typically,  the  smoke  from  individual  oii-weii  tires  combined 
and  rose  so  that  the  plume  top  reached  *'4.500  m.  often  with 
>hailow  convective  clouds  forming  on  the  windward  edge  of  the 
plume.  There  was  a  marked  vertical  wind  shear  leading  to 
differential  advection  of  the  plume,  with,  for  example,  the  lower 
pan  of  the  plume  being  transponed  southeasiwards  down  the 
Gulf,  with  the  upper  pan  moving  over  Iran. 

There  has  been  speculation  that,  because  of  radiative  self- 
heating,  aerosol  may  be  lofted'*  into  the  stratosphere,  where  its 
long  lifetime  would  allow  transpon  over  considerable  distances. 
Visual  observations  showed  the  plume  top  to  be  well  defined 
with  a  maximum  height  of  5.000  m.  On  one  or  two  occasions  a 
detached  smoke  layer  was  seen  a  few  tens  of  metres  above  the 
mam  top.  These  obserxations  in  springtime  show  that  the  smoke 
was  conrtned  to  the  lower  half  of  the  troposphere  and  significant 
amounts  were  not  being  lofted  into  the  stratosphere,  in  agree¬ 
ment  with  recent  assessments**' 

Particulates  and  chemistry  in  the  near-field 

Vertical  profiles  and  cross-wind  horizontal  runs  through  the 
plume  were  carried  out  on  28  March  1991.  120  km  downwind 
of  the  source  (dclined  as  29®  N.  48®  E)  south  of  Kuwait  city. 
This  distance  was  chosen  as  it  was  close  enough  to  the  source 
areas  that  all  emissions  could  be  measured  while  not  so  close 
that  the  inhomogeneity  of  the  individual  sources  rendered  the 
measurements  unrepresentative. 

-A  passive-caviiy  aerosol  spectrometer  probe  tPCASP)  was 
used  to  measure  smoke  panicle  concentrations  and  deduce  the 


1  Vertical  cross-secticn  of  the  tHume  made  in  a  aescenc 
mg  orohle  on  28  March  1991  l20km  from  the  source  a 
Smoke  density  a  SO}  tpoo.v  r.  c  Oj  fsoiiohnef  and 

NO.  idatfied  imei;  ct  0,  outsiOe  the  oKime  'sokd  >mel  and 
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table  I  Hyorocaroon  anatyses  of  plume  samples 

Samoie  A  Sample  6  Sample  C 

30  Marcr  191  30  Mwch  91  29  Marcn  91 


Mtimg 

rat« 

Einane 

Mifng 

ratio 

Ethane 

Mi««ne 

ratio 

Ethane 

Soec>«s 

IP  DO  V  * 

raiio 

IPPOv  1 

raiiO 

ippOw  > 

ratio 

emane 

25111 

1000 

5440 

1  000 

1  3X 

1000 

etneoe 

13  OX 

0  519 

1851 

0  340 

0034 

0026 

tKODane 

19952 

0794 

2.345 

0431 

0453 

0340 

prooenc 

1637 

0065 

0.220 

0  040 

0026 

0020 

Oenzene 

1057 

0042 

0.215 

0040 

0044 

0033 

toiue''e 

1  309 

0052 

0.X3 

0037 

0025 

0019 

5052 

0  201 

0409 

0075 

0013 

0010 

NMHC 

3X6 

1353 

4592 

8  44 

6  266 

471 

Analyses  of  1 6<iure  staintcss  steel  fiasns  Dy  gas  crwomaiograpriyfiame  Mxuiation 
oetec:ic^  Sample  A  most  poiiuteo  sample  120t(m  from  source  sample  B  typicai 
ooiiutec  sample  at  a  similar  Oistance.  sample  C  typical  tar-f«io  sample  (l.OOOkm 
f'om  sojrce 


mass  of  smoke  partiy-les  per  unit  volume  (smoke  density  !.  Con¬ 
tinuous  measurements  were  made  of  SO^,  NO»  and  ozone  (Oj) 
with  accuracies*  of  =5%,  and  =5%  respectively.  Stainless 
steel  flasks  were  pressurized  with  ambient  air  and  subsequently 
anaKsed  for  C;-Cj  hydrocarbons. 

Figure  I  shows  pollutant  concentrations,  together  with  the 
sum  of  0i  and  NO,  mixing  ratios,  from  a  descending  profile 
carried  out  alongwmd  through  the  centre  of  the  plume.  Also 
shown  IS  an  0%  profile  obtained  earlier  outside  the  plume.  The 
plume  had  a  complex  vertical  structure,  made  up  of  two  main 
layers  an  upper  one  with  a  top  at  4.600  m  and  base  at  ?.600m. 
and  a  deeper,  denser  lower  layer  with  a  top  at  3.200  m  and  base 
l.OOOm  The  lower  layer  was  capped  by  thin  (200mi  alto¬ 
cumulus  with  droplets  of  high  concentration  (270cm"*l  and 
small  effective  radius  (2.5  lami.  The  upper  ia.ver  was  adveaed 
eastwards  while  the  lower  layer  was  advected  southeastwards 
down  (he  Gulf,  trapped  b>  a  weak  temperature  inversion  at  the 
top  of  the  liquid  water  cloud 

Se\en  horizontal,  cross-wind  runs  over  the  same  ground  posi¬ 
tion  were  carried  out  at  1.420 m.  1,520m  (twice).  1.980m. 
r  440  m  and  4.420  m  (twicei  Figure  2a~d  shows  measurements 
from  a  run  through  a  .'0  km-wide  plume  (near  the  centre  of  the 
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fIG  2  ^y)rl7ontal  atKs-sectior  through  the  plume  at  2  aaO  m  on  28  March 
1991  a  Smoha  ttenstty.  0  total  scattering  coetfcient  et  0  48  pm  as 
measurea  By  an  mtagraiing  nephetomeier:  c  SO,.  0  0,  tsoM  Ime)  i»0. 
Short-aashea  Imel  ana  0,  .NO.  Itong-aasnea  ime' 
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FIG  3  Horizontal  run  at  an  attituoe  of  4  aoO  rn  through  the  area  of  snx)ke 
(-l.OOOkm  from  Kuwatt)  on  29  Marcr  1991  sryowing  a.  Oj  mixing  ratio.  0 
smoke  concentration 


lower  layer  at  2.440  m)  of  pollutants,  total  scattering  coefficient 
and  NO,  ■»‘Ot.  In  both  the  vertical  profiles  and  horizontal  runs 
all ihe  plume  indicators  (SO;.  NO,,  particulates  and  scattering 
coefficient!  were  highly  correlated;  this  correlation  was  observed 
over  several  days  of  measurements  close  to  the  source.  In  the 
lower  layer.  NO,  and  0%  mixing  ratios  were  negatively  corre¬ 
lated.  indicating  gas-phase  titration  of  O,  by  NO.  but  their  sum 
showed  little  horizontal  variation  (Fig.  2d),  indicating  little  loss 
of  0.11  either  by  reactions  with  smoke  panicles  or  in  the  oxidation 
of  SO;  to  sulphate.  In  the  upper  layer,  the  observed  correlation 
between  NO,  and  0>  was  zero  or  positive,  indicating  that 
photochemical  oxidant  production  had  already  been  initiated 

Analyses  of  hydrocarbon  samples  taken  in  different  pans  of 
the  plume  are  summarized  in  Table  1  The  flask  that  showed 
the  highest  mixing  ratios  of  hydrocarbons  (sample  A  i  was  filled 
at  an  altitude  of  2,600  m  in  the  densest  pan  of  a  smoke  plume 
on  30  March  1991,  120  km  from  source.  The  aggregate  of  the 
carbon  contained  in  all  the  identifiable  non-methane  hydro¬ 
carbons  (NMHC)  implies  the  presence  of  at  least  340p.p.b.\. 
of  carbon  as  unbumt  hydrocarbons.  Assuming  unidentified 
species  having  a  flame  ionization  response  are  hydrocarbons, 
we  infer  a  carbon  content  of  over  460p.p.b.v.  a  second  flask 
filled  on  the  same  day  (sample  B)  at  an  altitude  of  1.500  m  near 
the  bottom  of  the  plume  was  a  more  typical  sample;  compared 
with  a  UK  urban  air  sample,  the  maximum  mixing  ratios  of 
benzene  and  toluene  in  the  plume  sample  were  lower,  but  that 
of  n-hexane  was  substantially  greater. 

Two  supersaturation  cloud  condensation  nuclei  (CCN) 
spectra  taken  on  28  March  1991.  one  just  below  the  top  of  the 
upper  smoke  layer  at  4.500  m  and  the  other  in  clear  air  at  5.200  m. 
indicate  the  plume  contains  an  order  of  magnitude  more  CCN 
for  a  given  supersaturation.  For  example,  at  0.6%  super- 
saiuraiion.  CCN  concentration  increased  from  140  cm~''  to 
2.000  cm‘\  In  the  thicker  pans  of  the  plume.  CCN  were  so 
numerous  ( >4,000  cm~’)  that  the  counter  became  overloaded 

Far-AeM  measurements 

On  29  March  1991.  satellite  images  showed  smoke  over  central 
Saudi  Arabia,  some  distance  from  Kuwait,  oriented  southwest- 
northeast  and  being  advected  to  the  nonheast.  The  Oj  mixing 
ratios  and  panicle  concentrations  observed  in  a  horizontal  run 
through  this  area  of  smoke  are  shown  in  Fig.  3d  and  b.  In  this 
more  distant  plume,  it  is  likely  that  the  ozone,  enhanced  by 
40-60  p.p.b.v.  above  the  already  high  values  in  the  adjacent  clear 
air.  ts  formed  from  ozone  precursors!  NO,  and  NMHC)  released 
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from  the  oil  fires.  Sample  C  (Table  I),  taken  tn  this  distant 
plume,  has  a  ditlereni  distribution  of  hydrocarbons  from  the 
hydrocarbon  samples  collected  near  the  fires.  It  still  has  rela¬ 
tively  high  mixing  ratios  of  the  longer-lived  hydrocarbons,  such 
as  ethane  and  propane,  but  short-lived  hydrocarbons  such  as 
ethene  and  propene  are  substantially  reduced.  The  low  mixing 
ratios  of  shorter-lived  species  suggest  that  considerable  OH 
degradation  has  occurred,  as  is  consistent  with  the  enhanced 
ozone  mixing  ratio.  Although  the  NO,  mixing  ratios  in  the 
plume  were  below  the  limit  of  detection  (Sp.p.b.v.),  further 
ozone  production  remains  possible  because  mixing  ratios  of 
only  several  parts  per  10*^  (p.p.t.v.)  are  required  to  generate 
ozone*.  The  absence  of  many  of  the  more  reactive  and  therefore 
short-lived  hydrocarbons  means,  however,  that  further  ozone 
production  will  proceed  more  slowly. 

The  observations  of  this  mature  area  of  smoke  suggest  an 
interpretation  of  some  of  the  clear-air  observations  in  the  near 
field.  During  many  of  the  vertical  profiles  made  near  the  source 
region  but  outside  the  visible  plume,  ozone  mixing  ratios  much 
greater  than  the  normal  free  tropospheric  background  (30- 
40  p.p.b.v.)  were  encountered  (for  example.  Fig.  Id),  associated 
with  slightly  increased  panicle  density.  The  enhanced  ozone 
mixing  ratios  observed  in  the  clear  air  were  probably  formed 
from  NO,  and  NMHC  in  smoke  emitted  several  days  earlier 
that  remained  near,  or  returned  to.  the  source  area. 

Morphology  of  smoke  particles 

During  flights  in  the  smoke  plume,  ambient  air  was  drawn 
through  polycarbonate  filters  which  were  subsequently  studied 
with  a  scanning  electron  microscope  (SEM).  Figure  4a  shows 
a  typical  smoke  panicle  from  a  filter  exposed  180  km  from  the 
source,  on  30  March  1991.  The  smoke  panicles  arc  composed 
of  spherules  of  diameter  -0.1  um,  which  form  chains  of  hun¬ 
dreds  in  number  and  several  micrometres  in  length,  as  observed 
in  other  studies*  Panicles  on  filters  exposed  more  than 
1.000 km  from  Kuwait  (Fig.  46)  have  a  different  appearance 
from  those  from  the  near  field;  although  the  spherule  size  is 
similar  in  both  cases,  the  spherules  in  the  far  field  are  more 


TABLE  2  Estimates  Ql  SUength  ol  oil-well  sources 


Comoorcni 

This  stuov 
(Ml  vt''» 

Ref  2 
(Mt  yf‘‘) 

Ref  1 
(Ml  yr'M 

Ref.  3 
(Ml  yT"‘ 

Total  burn  rate 

202  5 

80 

63 

161 

Camon  as  gas 

161  2' 

60 

63 

Fine  carbon 

6  4-'. 

5 

58 

16.1 

particulate 

Sulphur 

61$ 

2 

_ 

Mitrogen  ifi  oxides 

0  42- 

05 

— 

— 

Estimates  oerwefl  from  ot>servaiions  ano.  for  companson.  those  assumed  t>y  recent 
assassmer^U^ The  mass  of  hyorogen  containad  m  HjS  arm  hydrocaroons  in  me  fuel 
has  been  taken  mto  account  wnen  aetermmtf^g  me  afferent  components  of  the 
emissions. 

*  Oeoosits  from  large  partially  bur r«o  od  droplets  near  the  well  heads  may  reduce 
this  mtai. 

t  Assuming  the  same  emission  factors  for  the  southern  and  northern  fees. 

t  This  quantity  coirid  be  up  to  2-3  times  higher  with  a  commensurate  reduction  m 
carbon  as  gas. 

f  Assuming  a  1.5%  sulphur  content  for  me  northern  fields  which  contribute  15% 
oueraN  to  the  emissions. 

likely  to  be  closely  packed  into  near-spherical  clusters.  Simitar 
behaviour  has  been  noted*  in  laboratory  experiments,  in  air 
aged  at  high  humidity. 

The  PCASP  determines  the  number  and  size  distribution  of 
particles  by  passing  them  through  a  laser  beam  and  measuring 
the  scattered  light  intensity,  using  a  relationship  between 
intensity  and  size  derived  from  experiments  with  latex  spheres 
of  known  radius.  The  smoke  density  is  calculated  by  integrating 
the  size  distribution  assuming  II)  that  collection  efficiency  is 
100%.  (2)  that  the  nonspherical  particles  in  the  plume  scatter 
the  same  amount  of  light  as  spheres  of  the  same  mass  (we 
conclude  from  previous  experimental  and  theoretical'^  work 
that,  for  the  particle  shapes  encountered,  this  assumption  is 
reasonable)  and  13)  that  the  density  of  the  smoke  panicles  is 
I  g  cm‘^.  in  line  with  Sokolik‘^  Suppon  for  this  approach  comes 
from  the  mass  loading  of  the  filter  samples,  which  indicate  that 
the  absolute  determination  of  smoke  density  from  the  PCASP 
is  not  grossly  in  error,  and  from  the  linear  relationship  between 


FKj.  4  Scanning  electron  m<rograohs  of  smoke  particles  coHected  on  fitters 
exposed  m  the  plume  on  a,  30  March  1991.  180  Km  from  the  source,  and 
0.  29  March  1991,  -1.000  km  from  the  source  fhe  filters  were  rrKAinted 
on  an  extendibie  boom  0  4m  from  the  arreraft  skm  and  the  air  sampled 
tsokineticaHy^^  The  dark  holes  m  the  photograph  are  the  pores  in  the  filter. 
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SO-  concenirauon>  and  smoke  densiu  isce  bclowi  which 
confirms  chat  the  instrumental  response  to  smoke  densitv  is 
linear  Nevertheless,  the  unceridiniies  in  these  assumptions  lead 
us  to  believe  that  the  absolute  values  of  smoke  densitv  could 
be  up  to  2-Mimes  as  high  as  those  shown  The  relative  variations 
tn  the  plots  are  stiH  valid 

Emission  estimates 

The  amount  of  sulphur  emitted  was  estimated  from  the  flux  of 
sulphur  through  a  crosswind  plane,  determined  from  the  product 
of  sulphur  concentration  and  the  wind  component  perpendicular 
tr  the  plane.  Sulphur  dioxide  measurements  from  the  horizontal 
runs  and  vertical  profiles  through  the  plume  were  used.  logelher 
w nh  the  winds  measured  bv  the  aircraft.  Above  5.000  m.  where 
the  SO:  monitor  will  not  operate.  PCASP  measurements  were 
used  to  derive  prow  SO:  values  through  the  linear  relationship 
e\  idem  in  Figs  I  and  2. 

On  28  March  1991.  when  virtuallv  all  the  plume  •'120km 
irom  the  southern  oilfields  was  sampled,  the  csiimaicd  sulphur 
emissions  were  5.“  Mt  yr'‘  Assuming  a  value  of  3.3%  for  the 
sulphur  content  of  the  southern  oilfields  (  R.  D,  Small,  personal 
communication),  this  source  strength  corresponds  to  an  oil 
burning  rate  of  r2.1  Mt  yr"‘.  We  assume,  based  on  previous 
e.vumates"  rha:  the  southern  oilfields  account  for  85%  of  the 
total  burning  rate;  this  implies  a  total  rate  of  202.5  Mtyr  ‘  or 
3  9  X  10’'  barrels  per  dav.  An  estimate  for  30  March  1991  gives 
a  similar  burning  rate.  Our  subiective  estimate  of  the  error  in 
this  source  strength,  based  on  the  temporal  evolution  and  spatial 
'  jrubilitv  over  the  h  measurement  period,  is  x40%. 

L  sing  the  relationships  between  concentrations  of  SO;.  NO. 
and  particulates,  we  estimate  the  source  strength  of  NO.  and 
namculates  from  the  sulphur  source  strength  tTable  2i  C  arbon 


•iG  5  iiiustratior  tne  «Mect  ot  the  smotse  on  oownwellmg  radiation  dur«n| 
a  cfosswtrvd  run  through  the  cemre  of  the  plume  a.  Smoke  density  b 
5nort«ave  fiun;  c.  infrared  fiu*  and  radiance  at  11  urn.  R..  Im  units  o' 
-ivv  rr.'‘  steradians''  a  the  zenith  microv.ave  onghtness  tem. 

De'stu'es  at  89  ano  157  GHz 


as  gas  ithe  mass  of  carbon  contained  in  CO.  CO-.  CHj  and 
N'MHC  •  is  estimated  from  the  residual  mass  after  the  mas.^ 
emissions  of  sulphur,  nitrogen  and  paniculate  carbon  are  sub¬ 
tracted  from  the  total  mass  emitted  and  an  allowance  is  made 
for  the  hydrogen  content. 

Effect  of  the  smoke  on  radiation 

Figure  5  shows  the  effect  of  the  smoke  on  the  radiation  during 
a  horizontal  run  on  31  March  1991  through  the  centre  of  the 
plume.  7tX)  m  below  its  top.  -“65  km  from  source.  The  total 
shonwave  (0.3-3  M-m)  dow'nwelling  flux  was  reduced  from  its 
clear-sky  value  of  800  Wm"'  outside  the  plume  to  zero  in  the 
centre,  and  corresponding  changes  were  seen  in  the  infrared 
and  microwave  regions. 

The  single  scattering  albedo.  a>„.  of  the  smoke  was  estimated 
in  two  independent  ws  -si.  we  measured  hemispheric 
reflectances  of  the  plu.  •  .u,.  using  upward-  and  downward- 
looking  pyranometers.  to  be  5  8%  over  a  range  of  solar  zenith 
angles  (25‘-48®).  The  Ginzburg  and  Sokolik*'  method  shows 
that  these  reffectancei>  correspond  to  a  value  of  u>(,  in  the  range 
0.5-0.55.  but  with  an  uncertainty  of  rO.l  because  of  assumptions 
made  about  the  scattering  properties  of  oil  smoke.  Second,  we 
measured  volume  absorption  coefficients  from  filters  exposed 
in  the  plume  using  the  integrating  sandwich  techniQuc''*.  and 
these,  when  combined  with  the  ncphelometer  scattering 
coefficients,  give  values  for  ti*,,  m  the  range  0.55-0.70(±0.l ). 
Thus  the  likely  range  of  u>,.  is  0. .50-0.65.  which  is  higher  than 
previously  published  values  for  oil  smoke"  but  lower  than 
that  of  0.8- 1.0  measured  in  non-carbonaccous  aerosols''. 

The  broadband  optical  depth  of  the  smoke  wa.s  estimated 
from  the  change  in  downwelling  shonwave  flux  as  the  aircraft 
ascended  or  descended  through  the  plume  Figure  6  show-s  the 
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-;G  6  Groao-aano  snor.-^ave  optical  oeoths  ot  tf'e  plume  *rom  •.r'-?e 
3ifferent  ascenamg  or  oescending  profiles  flown  through  the  o«ume  as  3 
'Lhction  3-  the  total  column  density  of  smoKe  m  the  oath  '-e  orot'ie  O'' 
29Marcn  1991  «as  further  away  from  '.he  source  ano  lower  smo-e  oensi:  es 
-vere  encountereo 

derived  optical  depth  as  a  function  of  smoke  column  densitv 
■'or  three  profiles  through  the  plume  top.  Solar  heating  rates  in 
;he  plume  were  calculated  from  shonwave  rtu.x  divergences  ;n 
.1  profile  at  local  noon  on  50  March  1991.  NO  km  from  the 
'ource  The  heating  rite  was  (>ptcallv  -50Kd“‘  throughout 
the  depth  f  5. .*00  m  i  of  the  smoke.  It  was  nearly  constant  :n  :he 
plume  because  the  smoke  density  increased  from  the  top  to  t.he 
.tase,  At  a  distance  of  550  km  from  the  source,  the  smoke  densitv 
vas  a  factor  of  four  lower  and  the  heating  rates  wrcre  50  K  d’ 
This  compares  wiih  typical  clear-sky  values  of  1  Kc'‘ 

Figure  5c  shows  the  downward  hemispheric  infrared 
i'O  um )  rtu’ies  to  be  95  Wm  '•  greater  in  the  centre  of  the  plume 
than  In  clear  air  and  the  vertical  downward  radiance  at  II  I  ti.rrt 
to  be  Increased  by  50  mW  m '  *  sieradians  ’’  *cm‘'  »■'  These  vvere 
measured  by  a  pyrgeometer‘*  and  a  narrow  held  of  view  1 1 
radiometer'"  respectively.  The  peak  in  radiance  at  95  km  along 
the  run  is  probably  due  to  cloud  above  and  is  also  seen  in  the 
microwave  measurements  (Fig,  5<I>  The  broadbano  infrared 
optical  depth  for  a  vertical  path  through  the  smoke  was  calcu¬ 
lated  in  a  similar  manner  to  the  visible  optical  depth.  It  was 
estimated  to  oc  0,34  =  0.05  for  the  profile  on  50  March  199| 
The  11.1-u.m  radiometer  was  not  operated  in  this  prorilc.  but 
the  maximum  vemcal  optical  depth  calculated  from  the  honzon- 
'al  run  in  Fig.  5c  was  0.30  =  0  02.  which  was  "5%  of  the  corre¬ 
sponding  broadband  optical  depth.  The  infrared  cooling  rate  :n 
the  plume  was  calculated  from  the  flu.x  divergences  to  be 
5  1  K  d'*.  which  is  1.0  K  d''  greater  than  that  in  adiaceni  clear 
Jir  Assuming  smoke  densities  based  on  the  PC.ASP  measure- 
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menis.  we  used  ihc  optical  depth  10  deduce  a  •. isible  extinction 
cocrticieni  of  “-13  m  g  which  agrees  broadly  with  oilier 
mea:kuremeni5  ot' oil  smoke  '  .  The  1 1  I  um  infrared  absorption 
coeiF.cient  was  deduced  to  be  0.7- 1 .3  m*  g  ' .  in  agreement  with 
laboratory  measurements*'* 

Figure  5J  shows  the  effect  the  plume  had  on  downward 
radiances  lor  brightness  lemperaiuresi  at  ^9  and  157  GHz, 
measured  by  a  microwave  radiometer.  .Atmospheric  emission  at 
these  frequencies  is  principally  controlled  by  water  vapour  and 
liquid  water.  The  increase  in  brightness  temperature  in  the  plume 
at  157  GHz  1 52  K I  is  three  limes  the  increase  at  .S9  GHz  1 S  K 1, 
out  both  are  well  correlated  with  the  infrared  radiance.  This 
correlation  vvas  seen  on  several  occasions  when  there  was  no 
trace  of  liquid  water  in  the  plume.  The  microwave  radiance 
increase  is  surprisingly  high  and  cannot  'oe  explained  bv  '.an- 
aeon  in  water  vapour  wuhin  the  plume  nor  bv  the  observed 
levels  ol  SO;  .Measurements***  of  the  absorption  coefficient  oi' 
diesel  soot  at  94  GHz  are  over  an  order  of  magnitude  too  small 
to  account  for  our  results. 

Discussion 

Our  best  estimate  of  oil  burning  rate  is  greater  than  those 
previously  assumed'  ’.  The  greater  CO;  emissions  that  this 
implies,  however,  represent  only  3%  of  total  annual  fossil-fuel 
emissions  and  would  still  have  a  negligible  etfeci  on  giobai 
climate  through  the  greenhouse  etfeci  The  sulphur  emission- 
calculated  arc  considerably  higher  than  those  taken  by  Brow  mng 
e(  d/.*.  and  althougn  we  would  expect  the  acidity  of  rainfall 
episodes  to  be  no  higher  than  previously  predicted,  because  it 
IS  limited  by  droplet  saturation,  these  episodes  will  extend 
funher  downwind  across  southern  .Asia.  As  predicted',  increased 
photochemical  oxidant  concentrations  were  observed  both  c:o>e 
-.0  the  >ource  and  further  afield.  But  the  emi^ilon^  of  \0  .ire 
smaller  than  those  previously  assumed’,  for  this  <eason  .iiu: 
because  they  are  rapidly  removed,  they  are  unlikely  to  influence 
the  global  .SO,  background  or  global  tropospheric  ozone 

The  estimated  smoke  emission  rate,  even  w  tth  our  barge  uncer¬ 
tainties.  is  similar  to  the  range  used  m  previous  assessments 
Furthermore,  the  product  of  our  derived  -isible  extinctior: 
coefficient  and  smoke  emission  rate  uhe  factor  that  determine 
the  radiative  impact  of  the  plume  for  a  given  meteorological 
sMuaiioni  IS  independent  of  the  uncertainties  in  our  estimate  of 
smoke  density  and  is  similar  to  that  inferred  :rom  Browning  c.' 
Or*  Previous  assessments*  assume  that  the  smoke  dia  not 
scatter  radiation,  the  scattering  that  we  observed  would  be 
insufficicni  10  alter  their  conclusions  significantiy.  The  observed 
infrared  effects,  although  previously  ignored.  wiU  only  very 
slightly  ameliorate  the  predicted  daytime  surface  cooling.  Our 
measurements  thus  confirm  the  assumptions  used  m  recent 
j>5esjmen(s'**\  and  the  observation  that  significant  amounts  of 
smoke  were  not  present  above  5.000  m  bears  out  their  predic¬ 
tions.  This  study  provides  strong  support  for  their  conclusions 
that,  although  the  ctfeas  may  be  significant  on  a  regional  scale, 
those  on  a  global  scale,  including  the  Asian  summer  monsoon, 
are  likely  to  be  insignificant.  — 
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THE  ESTIMATION  OF  THE  OPTIMAL  NUMBER  OF  PARTICLES  REQUIRED  FOR  A 
REGIONAL  MULTI- PARTICLE  LONG  RANGE  TRANSPORT  AND  DISPERSION  MODEL 


R.  H.  Maryon 

Meteorological  Office 
Bracknell, 

Berkshire,  UK 


1.  INTRODUCTION.  Simulation  of  the  evolution  of  the  Chernobyl  plume  over  the  10  or 
more  days  of  the  release  highlighted  the  influence  of  the  changing  meteorology  in  distorting 
the  plume  and  dispersing  materisd  over  extended  periods.  If  a  complex  plume  of  this  type 
is  to  be  simulated  by  a  Lagrangian  Monte  Carlo  technique — that  is,  by  the  release  of  very 
large  numbers  of  ‘particles’  into  the  model  atmosphere  at  source,  and  allowing  the  model 
winds  and  turbulence  parametrization  to  effect  the  spread — the  question  arises:  how  many 
particles  is  it  necessary  to  release  in  unit  time  to  simulate  the  plume  adequately?  It  is 
‘obvious’  that  the  more  particles  released,  the  greater  the  accuracy.  However,  even  with 
modem  computing  facilities  there  are  constraints  upon  the  number  of  particles  it  is  practicable 
to  utilise:  if  N  particles  are  released  every  hour  the  model  will  be  adverting  240iV  at  the  end 
of  10  days'  continuous  release,  with  potentially  very  high  storage  and  CPU  requirements.  On 
the  other  hand,  excessive  economies  in  the  number  of  particles  released  may  result  in  a  totally 
inadequate  plume  simulation.  What  is  the  optimal  configuration? 


Theory  can  help  a  little,  but  does  not  provide  a  clear  answer  in  situations  of  great 
complexity.  For  example,  a  standard  result  used  for  Lagrangian  particle  models  is  that,  in  a 
given  grid-cell. 


(1) 


where  c  is  the  ‘real’  concentration  in  a  grid-cell,  Q  is  the  mass  of  material  released,  V  is  the 
grid-cell  volume  (taken  as  constant  in  this  paper)  and  Vc  is  the  standard  deviation  of  the 
estimate.  Thus  the  ratio  A  is  a  measure  of  accuracy  (the  coefficient  of  variation),  which  is 
proportional  to  both  and  c~'/*;  for  a  given  value  of  A,  N  is  obtainable  from  c,  or  vice 

versa.  Ve/Q  in  (1)  can  be  interpreted  as  the  expected  fraction  of  particles  in  a  grid-cell  and 
the  complete  term  in  brackets  the  reciprocal  of  the  expected  number  of  particles.  From  (1) 
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confirming  riqsid  improvement  of  accuracy  with  particle  number  at  small  N  but  only  a  slow 
improvement  at  Urge  numbers. 


Expression  (1)  is  derived  from  elementary  probability  theory,  and  uses  the  initial  as¬ 
sumption  that  the  probability  of  a  particle  occurring  in  a  grid-cell  it  m/Q,  where  m  is  the 


‘real’  mass  occupying  the  grid-cell.  In  other  words,  (1)  supposes  that  the  diffusion  of  the 
plume  is  being  accurately  modelled.  ‘Real’  here,  is  taken  as  a  definitive  simulation  or  ‘valida¬ 
tion’  produced  by  a  Monte  Carlo  model  luing  an  indefinitely  large  number  of  particles.  Thus 
there  is  a  ‘correct’  number  of  particles  occupying  each  grid-cell,  which  may  be  expressed  as  a 
ratio  of  the  total  number  of  particles  released. 

Now  with  Monte  Carlo  models,  diffusion  is  effected  by  a  random  forcing,  or,  for  long 
range  transports,  quite  adequately  by  a  random  displacement  of  the  type 

*t+l  =  -I-  y/lLtKr  (2) 

where  xt  is  a  displacement  vector  at  time  t,  u  the  wind  velocity,  K  a  difiusivity  and  r  a 
random  number  generated  from  a  suitable  distribution.  Diffusion  will  not  necessarily  be 
modelled  consistently,  here,  for  extended  releases,  as  different  realizations  of  the  random  walk 
will  lead  to  relative,  local,  differences  in  particle  density,  which  will  be  acted  upon  by  the 
evolving  synoptic  wind  pattern.  Indeed,  different  realizations  of  the  diffusion  process  sire 
being  coupled  with  what  is  well  known  to  be  a  chaotic  system.  The  local  differences  will  be 
accentuated  with  the  passage  of  time.  Thus  simulations  of  a  complex  plume  using  relatively 
few  particles  will  lead  to  somewhat  different  configurations  from  the  hypothetical  validation, 
and  from  each  other.  The  object  of  this  exercise  is  to  estimate  the  number  of  particles  required 
to  remove  serious  error  due  to  these  effects.  Expression  (1)  tells  us  little  about  the  size  of 
N  required  to  model  a  complex  situation;  near  the  fringes  of  the  plume  c  will  be  indefinitely 
small,  and  with  Q  most  likely  unknown  it  is  very  difficult  to  arrive  at  a  sensible  estimate  for 
iV.  A  value  for  A  must  be  decided  upon,  and  even  a  choice  of  a  threshold  value  for  c  will  be 
species-dependent,  and  may  be  difficult  to  prescribe  a  priori  in  cases  of,  say,  nuclear  accident. 
If  c  and  Q  are  ignored,  and  attention  is  confined  to  the  spread  of  particles,  the  difficulty  arises 
of  defining  what  constitutes  an  acceptable  level  of  error. 

A  straightforward  way  of  estimating  the  number  of  particles  required  adequately  to 
simulate  a  complex  plume  of  many  days’  duration  is  to  carry  out  sensitivity  studies  leading 
to  an  empirical  result.  To  do  this  properly,  very  powerful  computing  facilities  are  needed, 
enabling  a  simulation  using  an  ‘indefinitely  large’  release  of  particles  to  be  used  as  validation. 
It  is  to  be  hoped  that  the  prior  assumption  of  the  adequacy  of  the  control,  for  a  period  at 
least,  after  release,  will  be  co:jfirmed  in  the  course  of  the  investigation.  Note  that  we  are  not, 
here,  trying  to  estimate  how  many  particles  are  needed  for  a  correct  answer  in  any  absolute 
sense,  as  no  model  is  capable  of  producing  such,  but  can  only  provide  an  approximation  which 
depends  upon  the  conditions,  the  model  configuration  and  parametrizations.  What  we  are 
trying  to  do  is  to  find  the  most  economical  configuration  which  produces  a  solution  which  is 
not  significantly  different  from  the  ‘best’  answer  that  the  multi- particle  model  is  capable  of 
producing. 

2.  THE  MODEL.  The  model  used  is  the  U.K.  nuclear  accident  response  model,  acronym 
NAME,  a  Lagrangian  3- dimensional  multi-level,  multiple  particle  model  of  which  a  brief  pre¬ 
liminary  description  is  contained  in  Maryon  and  Smith  (1989).  The  model  winds  and  other 
meteorology  are  obtained  from  the  U.K.  Met  Office  operational  ‘Limited  Area’  NWP  model. 
The  lowest  analysed  level  Is  the  atmospheric  boundary  layer,  which  is  derived  from  the  NWP 
profiles,  and  of  variable  depth.  The  upper  model  levels,  and  indeed  its  grid,  correspond  to 
that  of  the  underpinning  NWP  model.  The  resolution  is  0.4425  deg  both  lat  and  long,  but  the 
coordinates  are  rotated  so  that  Europe  benefits  from  the  relatively  uniform  mesh  of  Equatoriad 
regions. 

Particles  are  transported  by  the  model  winds,  including  the  analysed  vertical  compo¬ 
nent.  A  10m  wind  is  diagnosed  which  is  consistent  with  Monin-obukhov  similarity  theory.  All 
wind  components  are  zero  at  the  surface.  Sub-grid  turbulence  b  limulated  by  the  addition 
of  a  random  comptment  as  in  (2):  in  addition,  at  each  timestep  (15  min — roughly  the  eddy 


turnover  time  for  the  convective  boundary  layer)  all  the  particles  below  the  capping  inversion 
are  randomly  reassigned  in  the  vertical,  within  the  boundary  layer.  Horisontal  and  vertical 
diffusion  are  reduced  above  the  inversion.  The  meteorological  fields  from  the  NWP  model, 
and  the  rainfall  rate,  which  is  drawn  from  various  sources,  notably  radar,  are  archived  over  a 
10-day  period,  and  are  continually  updated.  The  model  hat  many  features  which  cannot  be 
described  in  the  space  available  here:  it  is  hoped  to  publish  a  fuller  description  shortly. 

3.  EXPEBJMENTAl  METHOD.  The  investigation  was  carried  out  using  the  NAME  model 
on  the  U.K.  Met  Office  Hitachi  EXIOO  mainframe  computer.  The  10-day  roll-over  archive  was 
‘frosen’  at  a  convenient  time  so  that  repeated  integrations  could  be  carried  out  using  identical 
meteorology.  This  archive  alone  required  nearly  300  Mbytes  of  storage. 

For  the  present  study  all  particles  were  released  from  a  reasonably  central  point  50  deg 
N  20  deg  E,  from  a  notional  source  in  the  boundary  layer,  the  release  commencing  12GMT 
23rd  June  1991.  Plots  of  the  particle  spread  show  that  over  the  10  days  some  drifted  out  of 
the  domain  and  are  lost  to  the  analysis.  Many  others,  of  course,  were  lifted  to  high  levels.  The 
particles  remaining  below  900  mb  were  used  in  the  calculations  for  this  study,  to  avoid  the 
processing  of  very  large  numbers  of  grid-cells,  and  to  concentrate  on  the  critical  near-surface 
particle  densities.  There  is  no  compelling  reason  to  assume  that  a  full  3-dimensional  analysis 
would  greatly  alter  the  results,  however. 

The  following  validation  runs  were  made: 

Exp  A:  4000  particles /hr  released  for  144  hr,  the  integration  continuing  up  to  240  hr. 

Exp  B:  as  Run  A,  but  6000  particles/hr  released. 

Exp  C:  12000  particles/hr  released  for  72  hr,  the  integration  continuing  up  to  240  hr. 

It  is  taken  on  trust,  for  the  moment,  that  some  or  all  of  these  runs  do  in  fact  constitute 
an  adequate  validation.  The  maximum  computer  core  used  was  38  Mbytes  and  the  longest 
integrations  extended  over  9  hr  CPU  time — despite  the  fact  that  in  the  interests  of  economy 
no  calculations  of  concentration  (as  opposed  to  particle  counts),  radioMtive  decay  or  wet 
and  dry  deposition  to  the  surface  were  made.  These  latter  diagnostics  we  all  dependent 
on  concentration,  and  if  calculated  may  be  expected  to  reflect  only  the  shortcomings  of  the 
particle  densities,  which  are  to  be  analysed  here. 

Identical  meteorological  data  were  employed  for  trial  integrations  corresponding  to  Exp 
A,  B  and  C  using  releases  of  2,4,8,16,32,64,128,256,512,  1000,2000,3000...  particles/hr,  cul¬ 
minating  in  a  release  of  similar  size  to  the  validation  but  using  a  different  realisation  of  the 
random  walk  (2).  It  was  not  deemed  necessary  to  include  all  the  larger  integrations  for  Exp 
C,  however. 

4.  ANALYSIS  OF  RESULTS.  A  suitable  statistic  for  representing  the  error  of  the  trial  runs 

^  =  (») 

where  t  sums  over  all  bottom  layer  grid  cells,  nia>  are  the  numbers  of  particles  in  grid-cell 
t  for  releases  of  Wa,  IVs  particles/hr,  over  a  period  of  t  hr,  where  subscripts  a,  b  indicate  the 
validation  and  trial  runs  respectively. 

A  digrcMioni  the  statistic.  It  was  noted  that  stsdistic  S  is  not  far  removed 

from  a  x-Ml<tare  statistic,  if  it  is  assumed  that  the  grid-cell  contents  of  the  validating  runs 
constitute  an  ‘expected  distribution’  of  the  particles.  For  air  concenttatum  a  suitable  statistic 
can  be  defined 

x2 


with  the  RHS  subscripts  as  in  (3).  In  terms  of  the  particle  analysis  this  becomes 


X?  = 

if  V  is  constant.  Here,  R  is  the  ratio  It  it  not  practicable  to  use  this  formula,  which 

it  a  function  of  Q.  However  the  summation 

E 

i 

is  a  closely  related  statistic  where  Rna  is  the  weighted  nnmber  of  particles  in  the  trial 
grid-cells — we  are  saying,  in  effect,  that  each  trial  particle  corresponds  to  R  of  the  validation 
particles,  which  in  the  trial  runt  are  forced  to  ‘adhere’  due  to  the  loss  of  degrees  of  freedom 
where  fewer  particles  are  avaulable.  Thus  in  (4)  we  compute  the  squared  difference  between 
this  weighted  nnmber  of  trial  particles  and  the  number  of  validation  particles  in  the  grid-cells. 

A  test  might  now  be  east  entirely  in  terms  of  particle  numbers,  without  dependence 
on  unknowns  such  as  Q.  The  degrees  of  freedom,  v,  of  the  test  is  taken  as  the  number 
of  grid-cells  occupied  by  S  or  more  particles  of  the  validating  plume;  those  with  4  or  less 
are  combined  in  one  ‘bin’.  One  possible  level  of  significance  for  the  test  is  50%,  for  which,  at 
random,  half  the  acceptable  realisations  would  have  Xp  >  Xso-  degrees  of  freedom  for  such 

a  test  greatly  exceed  the  values  usually  published  in  table  form,  but  to  a  close  approximation, 

(see  Dixon  and  Massey,  1969;  for  other  significance  levels  the  formula  has  an  additional  term 
requiring  the  normal  deviate).  Clearly,  as  s'  co,  Xso  “*  *'•  This  limit  applied  to  (4)  leads 
to  the  interesting  deduction  that  in  individual  cells  t  we  are  looking  for 

ttio  “  ^  Rfti^  ^  71, 'q  -f-  (6) 

(to  a  rough  approximation),  or  transposed  for  air  concentration 

’*io  ^  /tSia  V  Tljp  /n^a  V 
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for  a  successful  trial. 

Are  the  validations  inadequate?  Fig  1  is  a  plot  on  logarithmic  axes  of  against  the 
number  of  particles  released/hr  for  Exp  B.  The  three  curves  correspond  to  the  situation  at 
72,  144  and  240  hr  into  the  integration;  the  continuous  release,  it  will  be  recalled,  extended 
to  144  hr.  The  approximation  of  the  curves  to  straight  lines  of  unit  slope  is  noteworthy.  The 
plots  imply  oc  l/N — interpreting  (4)  as  a  form  of  squared  error  this  it  a  clear  manifestation 
of  the  N~^l^  re^me  discussed  earlier. 

Where  this  relation  hold*  the  error  due  to  reetricted  particle  nambers  is  predominant. 

However,  as  the  particle  numbers  increase,  the  slopes,  at  first  imperceptibly,  and  towards 
the  extremity  more  noticeably,  depart  fimn  unity  (the  Ml  hand  extremes,  where  the  plume  is 
ofttdely  discretised  by  a  handful  of  particles,  need  not  concern  us).  Where  very  large  nmnbers 
of  partidee  am  employed,  the  deoease  In  the  reduction  of  error  with  nmnber  refleeU  the 
inflaettoe  of  afferent  teaHemtione  of  the  d^feeion  process  (2).  Whichever  validation  we  look 
at,  and  however  many  particles  am  utUiaed,  a  departum  fimn  unit  slope  occurs  fm  this  reason. 
This  decrease  in  errm  redncti«m  prevents  any  does  ^proadi  of  x?  to  xio  **>7  other  of  the 
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conunonly  used  significance  levels)  even  where  a  trial  uses  as  many  particles  as  the  control 
run  (but  in  a  different  realisation).  LogXu)  for  the  three  curves  lies  between  2.48  and  3.49: 
if  the  unit  slope  sections  of  the  curve  are  produced  to  intersect  these  values  the  ordinate 
corresponds,  roughly,  to  the  numbers  of  particles  released  in  the  validation.  This  is  no  doubt 
due  to  the  fact  that  Xjo  fite  number  of  particles  released  are  both  closely  related  to  the 
number  of  grid-cells  covered  by  the  simulated  plume. 

The  implication  here  is  that,  in  these  validations  and  trials,  the  more  particles  released 
the  more  grid-cells  were  covered,  which  leads  us  to  question  the  adequacy  of  the  validations. 
Indeed,  plots  of  the  proportion  of  validation  grid-cells  covered  against  the  number  of  particles 
released  in  the  trials  (e.g.,  Fig.2  for  Exp  C)  fiul  to  show  (at  least  at  144  and  240  hr)  the 
levelling  off  for  large  numbers  of  particles  which  would  be  expected  if  the  validating  runs  were 
adequate,  and  the  trials  with  many  particles  reasonable  approximations  thereto.  Again,  plots 
(not  shown  here)  of  mean  proportionate  error  defined 

/jffnjs  — ni^\ 


(where  the  angle  brsuricets  imply  an  average  over  the  ensemble  of  grid-cells)  likewise  fail  to 
asymptote  sufficiently  as  the  validating  particle  numbers  are  approached,  except  perhaps  at 
72  hr,  and  this  only  close  to  the  validation. 

At  first  sight  these  seem  discouraging  results.  However,  where  a  random  walk  is  applied 
to  many  particles,  some,  after  a  number  of  timesteps,  will  be  displaced  a  considerable  distance 
&om  the  mainstream.  The  ensemble  of  accumulated  displacements  will  take  the  form  of  a 
Gaussian  distribution,  and  where  extremely  large  numbers  of  particles  are  deployed,  some 
displacements  will,  due  to  the  accumulation  of  random  effects,  spread  further  and  further 
into  the  tails  of  the  distribution.  The  more  extreme  displacements  may  eventually  come  to 
represent  a  virtually  insignificant  proportion  of  the  ensemble,  but  will  be  reflected  in  Fig.2 
because  some  of  them  will  spread  into  fresh  grid-cells,  and  in  the  proportionate  error  because 
grid-ceUs  containing  very  few  particles  will  carry  a  disproportionate  weight  in  expressions  such 
as  (7)  (via  its  denominator).  A  similar  argument  applies  to  expression  (4),  the  statistic. 

Now  with  extremely  large  numbers  of  partides  a  simulated  plume  can  expand  almost 
indefinitely;  thus  in  (2),  with  a  timestep  At  of  16  min  and  a  diffhsivity  JT  set  at  6000Tn^/s 
(if  anything  set  slightly  on  the  large  side  for  this  project  so  as  to  avoid  underestimating 
minimum  particle  requirements)  those  particles  repeatedly  being  displaced  by  near-maximum 
values  of  r  will  move  about  3km  from  the  mainstream  at  every  timestep.  As  a  fraction  of  the 
total  numbers  released,  the  number  of  particles  reaching  such  peripheral  grid-cells  would  be 
negligible.  The  tails  of  the  Gaussian  distribution  in  typical  short  range  dispersion  simulations 
are  taken  as  representing  very  uncommon  effects  of  the  turbulence  or  meandering  which 
may  occur  in  different  realisations.  Fw  a  tpecific  long  range  transport  sitaotson,  analogous 
extremes  will  occur  in  ixmumerable  places,  but  do  they  occur  aU  around  the  plume  at  all 
times  and  in  all  directions?  It  is  likely  that  parametrisations  such  as  (2)  are  not  everywhere 
realistic  when  applied  to  a  single  realisation,  or  at  least  that  in  reality  the  spread  of  material 
is  totally  insignificant  in  many  of  these  grid-ceUs.  The  fixed  difiusivity  used  here  will  itself 
be  unhelpful,  while  the  release  of  very  large  numbers  of  particles  may  even  exacerbate  the 
situation.  Parametrisittg  ff  as  a  function  of  the  meteorology  or  conditions  may  he^,  but  would 
not  solve  the  problem.  It  is  concluded,  then,  that  diagnostics  such  as  Fig.2  and  expression 
(7)  should  not  unquestioningly  be  taken  as  approj^iate  for  the  problem  in  hand. 

Soma  ampiiicsd  anppcni.  It  was  thus  necessary  to  adopt  a  different  ^prooch,  returning  to 
the  point  at  depart  ure,  the  error  statistic  5,  defined  in  (3).  Figs.3a,  3b  show  plots  at  S  against 
log  of  porticia  number  fiar  the  experiments  A  and  B.  It  will  be  obaarved  that  dUbrances  are 
very  sB^  iq»  to  ns  s  4900 — that  is,  error  for  a  vaildorion  of  4000  partklas/hr  is  very  dose 


indeed  to  that  for  6000,  even  after  240  hr  simulation.  The  largest  errors  in  Fig.3  are  attached 
to  the  72  hr  integrations,  no  doubt  because  of  the  relatively  large  values  of  n,'a  and  tin  which 
can  occur  at  relatively  small  times.  Normsdising  S  by  reference  to  the  validating  runs  defines 
a  statistic  T; 


(8) 


plotted  in  Fig.4  for  Elxp  B.  It  is  apparent  that  this  normalised  error  increases  with  time  into 
integration,  especially  for  transports  following  the  cessation  of  the  release,  when  the  statistics 
no  longer  benefit  from  the  near-source  plume.  Differences  between  plots  for  runs  A  and  B  are 
virtually  indiscernible  up  to  =  4000,  for  this  statistic. 


The  diagnostic 

Rn  = 

is  also  of  interest  in  that  the  ratio  for  a  perfect  simulation  would  be  imity,  and  it  demonstrates 
the  rate  at  which  the  trials  approach  a  good  representation  of  the  characteristic  scatter.  This 
is  illustrated  in  Fig.5  for  run  C.  The  curves  for  runs  A  and  B  converge  marginally  more 
quiclcly,  if  anything,  and  again  are  near-identical  up  to  4000  particles/hr.  It  is  necessary  to 
take  care  in  interpreting  some  of  these  statistics:  the  rate  at  which  these  curves  converge  gives 
an  erroneous  impression  of  the  number  of  particles  necessary  to  give  an  adequate  simulation. 

The  near  identity  of  runs  A  and  B  for  statistics  R„  and  T  and  (most  importantly)  the 
strong  similarities  of  these  runs  for  the  unmodified  error  5  imply  that  6000  particles /hr  should 
constitute  an  adequate  validation  for  the  sittiation  modelled,  contrary  to  the  impression  given 
by  Figs.2  and  3  where  the  results  are  influenced  by  fringe  effects  of  questionable  significance. 

To  check  on  the  significance  of  these  peripheral  effects  a  count  was  made  of  the  average 
number  of  particles  in  those  validation  grid-cells  which  the  trials  failed  to  resuh.  These  values 
are  plotted  in  Fig.6  for  run  C,  for  example,  with  the  control  release  of  12000  particles/hr 
(for  72  hr).  This,  for  the  most  part,  underlines  the  minute  number  of  particles  contained  by 
these  peripheral  cells;  a  cell  with  ten  particles  holds  less  than  0.1%  of  an  hour’s  release,  in  a 
region  which  might  well  not  be  affected  in  reality.  Of  course,  a  minority  of  cells  will  contain 
considerably  more  than  the  average:  for  example,  modelling  the  240  hr  plume  of  run  C  with 
512  particles  left  no  lest  than  60%  of  the  validation  cells  unreached  (Fig.2);  of  these  10  cells 
contained  40  -  50  particles,  8  cells  51  -  61  (the  nutximnm),  but  the  overall  average  was  5.86 
particles /cell.  What  proportion  of  this  very  wide  diffusion  was  realistic? 
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A  further  point  to  note,  here,  is  that  the  more  trial  particles  released,  the  more  grid- 
cells  are  reached  which  were  not  occupied  in  the  validation  (although  they  would  have  been 
in  a  hjrpothetical  simulation  with  an  indefinitely  large  number  of  particles!).  These  effects  are 
relatively  insignificant. 


5.  ESTIMATING  OPTDdAL  PARTICLE  NUMBERS.  It  remains  to  make  some  estimate  of 
the  optimal  release  in  the  given  set  of  circumstances,  although  it  must  be  emphasised  that 
there  is  no  perfect  or  absolute  answer  to  the  problem.  Two  further  diagnostics  were  calculated: 


1  ^  Stt  ~  ^n^^l 

5  dn  5n.t.iAn 


(10) 


and 
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In  these  expresaimu  the  5  terms  on  the  RES  are  as  in  (3)  with  the  subscripts  reftning  to 
consecutive  trisUs  (tee  Section  3)  or,  for  tub(9),  the  final  trial,  which  uses  a  similar  number 


of  particles  to  the  validation  but  with  a  dilFerent  realization  of  the  random  walk.  Thus  S\ 
represents  the  proportionate  decrease  of  error  with  particle  number  between  consecutive  trials, 

5]  the  proportionate  decrease  of  error  between  individual  trials  and  that  of  the  eijui-particle 
realization  (an  ‘irreducible  error’).  5 1  is  plotted  in  Fig.7  against  the  mean  of  n  and  n  + 1  (and 
in  Fig.8  St  against  the  mean  of  n  and  q),  and  for  plotting  purposes  expressed  as  percentage 
improvement  per  release  of  10  particles/hr. 

These  plots  are  notable  for  the  degree  to  which  the  curves  for  the  different  periods 
of  integration  coalesce.  Fig.7  suggests  that  the  improvements  obtidnable  afta  the  512/1000 
particle  step  are  very  slight.  The  strong  vertical  in  Fig.8  is  to  be  expected,  as  the  abscissa 
changes  little  for  triab  releasing  few  particies.  Nonetheless,  this  somewhat  unorthodox  ab¬ 
scissa  highlights,  at  the  particularly  sharp  bend  in  the  curve,  the  stage  at  which  increasing 
the  particle  numbers  produces  decreasing  dividends  in  accuracy.  The  implication  of  all  the 
diagrams  from  Fig.3  onwards,  and  of  a  number  of  other  diagrams  and  statistics  for  which 
there  is  no  room  here,  is  that,  given  that  a  reasonable  validation  for  integrations  of  72  up 
to  240  hr  was  obtained,  the  rate  of  improvement  is  rather  small  for  increases  beyond  500 
particles/hr,  although  Fig.5  suggests  that  1000  particles/hr,  if  feasible,  is  distinctly  prefer¬ 
able.  The  rate  of  improvement  is  very  small  after  1000  particles/hr.  Of  course,  the  very 
nature  of  the  plume  changes  between  72  and  240  hr,  and  these  effects  are  not  sdways  apparent 
in  the  statistics:  attainable  accuracy  reduces  with  time.  If  improving  computer  capabilities 
enable  the  deployment  of  larger  numbers  then  so  much  the  better,  but  the  present  exerdse 
suggests  that  improvements  in  accuracy  wiU  be  marginal,  and  eventually  constrsuned  by  the 
workings  of  the  diffusion  parametrization.  The  use  of  a  fixed  (and  rather  large)  diffusivity  in 
this  investigation  was  not,  with  hindsight,  helpful. 

Choosing  a  threshold  value  for  particle  number,  below  which  a  grid-cell  is  ignored  for  the 
purposes  of  the  investigation,  would  lead  to  an  interesting  reappraisal  of  the  test  discussed 
earlier.  This  would  again,  however,  require  an  a  priori  dedsion  as  to  realistic  spread,  and  is 
not  pursued  here. 

Obviously  the  statistics  derived  here  are  situation-specific,  depending  on  the  period  of 
the  release,  and  particularly  the  presence  or  absence  of  denser  concentrations  around  the 
source.  The  meteorological  conditions  may  be  important,  particularly  if  the  situation  is 
changing  rapidly,  or  there  are  strong  vertical  velodties  involved.  Reruns  of  this  study  with 
other  model  configurations,  release  scenarios  and  meteorology  would  naturally  be  desirable, 
but  it  is  felt  may  not  change  the  condusions  significantly.  A  useful  mode  of  operation,  in 
practice,  may  be  to  split  partides  after  a  certain  lapse  of  time,  utilising  the  computer  storage 
of  partides  which  have  drifted  out  of  the  domain. 
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FIGURE  CAPTIONS. 

Flg.l.  E:q>.B:  A  plot  <m  logarithmic  axes  of  against  the  number  of  partides  releaied  per 
hour  Cor  the  trials.  The  72,  144  and  240  hr  integratkos  are  distingnished  by  drcles,  crosses 
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and  dots  respectively.  This  applies  to  all  the  diagrams,  which  also  all  have  identical  abscissae. 
Fig.2.  Exp.C:  percentage  of  validation  grid-cells  reached  by  the  trial  integrations  plotted 
against  the  log  of  the  number  of  particles  released/hr. 

Fig.3.  Error  statistic  5  as  defined  in  (3)  plotted  against  log  of  particle  number:  (a)  Exp. A 
(4000  particle/hr  validation),  (b)  (6000  particle/hr  validation). 

Fig. 4.  Exp.B:  Normalised  error  as  defined  in  (8)  plotted  against  log  of  particle  number. 
Fig.5.  Exp.C:  Ratio  of  RMS  errors,  expression  (9),  plotted  against  log  of  psuticle  number. 
Fig.6.  Exp.C:  average  number  of  particles  in  validation  grid-cells  not  reached  by  trials,  plotted 
against  log  of  particle  number. 

Fig.7.  Exp.B:  Proportionate  decrease  in  error  with  particie  number,  5i,  defined  as  in  (10), 
plotted  against  log  of  mean  of  n  and  n  -r  1. 

Fig.8.  Exp.B:  Proportionate  decrease  of  error  Sj,  defined  as  in  (11),  plotted  against  log  of 
mean  of  n  and  q, 
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Model  Intercomoarison  of  the  Kuwait  Smoke  Plume 
A  meeting  held  at  the  Meteorological  Office,  Bracknell  S*^  November  1991 
R.  Maryon’,  H.  ApSimon^  and  I.M.  Lowles^ 

1  The  Meteorological  0£Sce.  London  Road.  Bracknell.  Berkshire.  RG12  2SZ.  UK 

2  Ingterial  College  of  Science.  Technology  and  Medicine,  London.  SW7  2AZ.  UK 

Following  the  ARATS  2  meeting  (16-17  September  1991)  Imperial  College  and  the 
Meteorological  Office  have  agreed  to  study  the  dispersion  and  deposition  of  the  Kuwaiti 
smoke  plume  using  two  different  Lagrangian  multipaiticle  models,  namely  3-DRAW 
(Imperial  College)  and  NAME  (Met.  Office).  Both  organisations  have  agreed  on  the 
following  model  input  data: 

1.  Source  Positions: 

a) .  29.25  N  47.40  E  (derived  from  two  sources  at  29.3  N  47.6  E  and  29.2  N  47.2  E) 

b) .  30.50  N- 47.80  E 

2.  Source  Terms:  ^ 

The  source  located  at  position  a),  released  3.95  x  10*  tons  of  Carbon  per  year 
The  source  located  at  position  b).  released  6.58  x  10*  tons  of  Carbon  per  year 

3.  Release  Heights: 

There  are  to  be  two  release  heights,  one  between  1000-2000  m  (900-800  mb)  representing 
an  elevated  release  and  a  second  between  ground-level  and  1000  m  (10(K)-9(X)  mb) 
corresponding  to  a  release  within  the  boundary  layer. 

4.  Release  Details: 

The  release  will  start  at  00  GMT  27/3/91  and  will  continue  for  48  hours. 

The  models  will  release  5(X)  particles  per  hour  over  the  same  48  hour  period  (00  GMT 
27/3/91-00  GMT  29/3/91). 

5.  Meteorological  Data: 

Both  3-DRAW  and  NAME  utilise  meteorological  fields  derived  from  the  Met.  Office 
numerical  weadier  prediction  model  FINE  MESH.  This  model  covers  the  area  1 1.25  W  to 
109.6875  E  and  from  4.5  N  to  54.0  N  with  129  x  67  grid  points  at  intervals  of  0.9375 
degrees  longitude  and  0.75  degrees  latitude.  The  horizontal  components  of  the  windfield 
are  calculated  on  a  wind  grid  offset  half  a  cell  from  that  described  above.  However, 
vertical  wind  velocities,  temperature  profiles,  cloud  cover  and  rainfall  data  are  all 
calculated  on  the  nradel  grid. 


FINE  MESH  model  data  will  be  supplied,  by  the  Met.  Office,  for  the  period  00 
GMT  27/3/91  to  00  GMT  29/3/91. 

The  models  will  compare  the  following  results  for  1200  GMT  28-29/3/91. 

j.  Trajectory  end-points. 

a.  Ground  and  elevated  level  concentrations. 

iii.  Deposition. 

Vertical  concentration  profiles  of  the  smoke  at  29  N  48  E  on  the  28/3/91  will  be  calculated 
and  compared  with  airborne  measurements  performed  by  the  Met.  Office  flight. 

7.  Timetable: 

Meteorological  data  supplied  by  the  Met.  Office  -  end  of  November  1991 
Initial  model  runs  -  end  of  January  1992 
Revision  and  fine  tuning  -  end  of  February  1992 
Report  -  March  1992. 
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